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The work described in this report was performed by the Quality 
Assurance and Reliability Division of the Jet Propulsion Laboratory. 
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ABSTRACT 


This technical memorandum presents the computer program listing 
for the reliability block diagram computation program described in 
Technical Report 32-1543, Reliability Computation From Reliability Block 
Diagrams , Jet Propulsion Laboratory, Pasadena, Calif., Dec. 1, 1971. 
The program is written in FORTRAN V and is currently running on a 
UNIVAC 1108. Each subroutine contains a description of its function. 
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I. INTRODUCTION 


This technical memorandum presents the computer program listing 
for the reliability block diagram computation program described in 
Technical Report 3Z-1543, Reliability Computation From Reliability Block 
Diagrams , Jet Propulsion Laboratory, Pasadena, Calif., Dec. 1, 1971. 
The program is written in FORTRAN V and is currently running on a 
UNIVAC 1108. Each subroutine contains a description of its function. 
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ALL INTEGER FORMATS RIGHT JUSTIFIED. COLUMN NUMBERS APPEAR INSIDEMAIN 

13C 

15 
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PARENTHESIS ( ). ANY SERIES OF BLOCK NUMBERS STARTS AT THE LEFT 

MAIN 

140 

18 

c 

MOST COLUMN AN FILLS IN SUCCESSIVLLY TO THE RIGHT IN 12 FORMAT. 

MAIN 

150 
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MAIN 

160 
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c 


MAIN 

17 C 

19 
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BLOCK NUMBER(i-2)f IiNPUT TO THAT 3L0CK NUMBER ( 3- 32 ) t OUTPUT FROM 

MAIN 
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THAT BLOCK N UM& E R ( 33 -63 > . 

MAIN 

19C 
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MAIN 
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22 
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CARD IN THIS SERIES. 

MAIN 

210 

23 
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MAIN 

220 

24 
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NUMBER OF SENSE cL0CKS(l-21. 

MAIN 

230 

25 

c 

IF THERE are SENSE BLOCKSt SENSE BLOCK NUMBER(l-2)t STANDBY BLOCKSMAIN 

240 

26 

c 

CONTROLLED BY THAT SENSE EL0CK(3“6C). 

MAIN 

25C 

27 
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MAIN 

2 60 

28 

c 

NUMBER OF EQUIVALENT BLOCK SETSU-2). 

MAIN 

27C 

29 

c 

EQUIVALENT BLOCKS HAY NOT BE CONTROLLED BY A SENSE BLOCK « SW IT C H ) . 

MAIN 

2 80 

3C 

c 

IF THERE ARE EQUIVALENT BLOCKS* EQUIVALENT BLOCK NUMBERS < 1-4 C ) FCRMAIN 

290 

31 

c 

EACH SET. 

MAIN 

3 CO 

52 
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MAIN 

310 

33 

c 

DISTRIBUTION TYPE(l-2)f MISSION TIM£(3-14l E FORMAT. - 

MAIN 

32C 

34 

c 

1 r EXPONENTIAL# 2 NC T AVAILABLE. 

MAIN 

330 

35 
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MAIN 

34C 

36 

c 

ACTIVE PARAMETERS FOR EACH BLOCK OTHER THAN SENSE BLOCKS « SWi TCHE S 

)MAIN 

35 G 

37 

c 

RC IS SET EQUAL TO 1.0 IF(15-25) LEFT BLANK. 

MAIN 

3 60 

38 

c 

BLOCK NUMS£Ra-2)f LAMfcCA(3-14) E FORMAT# RC(l5-25> F FORMAT. 

MAIN 

37C 

33 

c 

LAST CARD IN THIS SERIES HAS LASTOGJ SET EQUAL TO 5. 

MAIN 

380 

4C 
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MAIN 

39 C 

41 

c 

IF NO DORMANCY INVOLVED# A BLANK CARD# OTHERWISE - 

MAIN 

4 CO 

42 

c 

3LANK«l“2>f dormancy FACT0K(3*14) E FORMAT. TC BE MULTIPLIED BY 

MAIN 

41C 

43 

c 

ACTIVE LAMBDAS OR 99(1-2)# 3LANK(3-8G) AND READ INDIVIDUAL DORMANTMAIN 

420 

44 

c 

LAMBDAS - BLOCK NUMBER(i-2)# LAMBDA DORM ANT ( 3-14 ) E FORMAT. 

MAIN 

43 0 

45 

c 

LAST CARO IN THIS SERIES HAS LAST(SG) SET EQUAL TO 6. 

MAIN 

440 

46 
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MAIN 

450 

47 

c 


MAIN 

4 60 

48 

c 

SWITCHING OPTIONS FOR EACH SENSE BLOCK. 

MAIN 

47C 

49 
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MAIN 

4 80 

5C 

c 

G NO SWITCHING( PRCbABlLlTY Cr SWITCH WORKING = l.C). 

MAIN 

49 

51 

c 

SENSE 3LCCK(i-2l. 

MAIN 

5C0 

52 
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MAIN 

510 

53 

c 

1 - CONSTANT PROBABILITY THAT SWITCH WORKS. 

MAIN 

5 20 

54 

c 

SENSE BLOCK 1 1-2 ) #BLANK( 3-14) #PRCBABILITY (15-25) F FORMAT# K8C ). 

MAIN 

53 0 

55 
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MAIN 

540 

56 

c 

2 - DORMANT FAILURE RATE FOR SWITCH. 

MAIN 

55C 

57 

c 

SENSE 3L0CKC1-2) #LAM3GA OORMA NT ( 3-14 ) E FORMAT# 2(30). 

MAIN 

5 60 

58 
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MAIN 

570 

59 

r 

3 - dormant and ACTIVE FAILURE RATE FOR SWiTCHi 2 C ARC S/SW I TCH ) . 

MAIN 

5 80 

60 
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SENSE GLCCK(l-2 ) flAMBDA DORM AN T ( 3 -14 ) E FORMAT# 3(80). 

MA IN 

590 

61 
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SENSE 3L0CK(l-2 J #LAMOOA AC T IV E ( 3-14 ) . 

MAIN 

600 

62 
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MAIN 

610 

63 
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4 = NOT available. 

MAIN 

6 2Q 

64 

c 

5 - NOT available. 

MAIN 

63C 

65 

c 


MAIN 

640 

66 

c 

LAST CARD IN DATA DECK FAS LAST(8C) SET EQUAL TO 7* 8# OR 3 - 

MAIN 

650 

S7 

c 

7 = RECALCULATE WITH NEW PARAMETERS# 3 - NEW DIAGRAM# 9 = END. 

MAIN 

6 60 

68 
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MAIN 

670 

69 

c ♦♦ 
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6 30 

70 
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MAIN 

69C 
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COMMON/ ALLSUB /M If M2 #M 3 #M 4# M 5# M8#M7#Ll#L3fL4#L5fL6 

MAIN 

7C0 

72 


CGMMGN/L6SUB/IB ( 15*2 #SC ) # IIS (50 ) * IRB (50* 30 )* IT BMP (41 #20 ) *N1# 

MAIN 

71G 

73 


1 IS AVE( 50 f 200 ) t ISUC# JSUC * JSl# JSA VE ( 5C # 20C ) . IJSI# TOT #PRC3( 8 5 )# 

MAIN 

7 20 

74 


E F RO ci ( 65 ) *IFRINT fUSTDB Y( 15 ) *NSBYMX * NGPRNT 

MAIN 

730 

75 


C0MMCN/SMSU3/NPUP ( 200 »N PUP MX # NS OUT (50) * NB OT MX * NB I N ( 50 ) # NB INMX • 

MAIN 

740 
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76 


1 NBN'JM(5C » »N6:NM AX 

MAIN 750 

77 


OIKENSIONi RC(50)*r4P(5C)#RT(65flGl)f PSNS (ICl) t IS 1 5C > *TL ( 53 1 

MAIN 76C 

78 


dime: NSION IN <15 ) tIGUT ( 16 ) f JR & (3C) fNRB (3C ) #KR B (30 tKS AVB( 5C#2CC ) 

MAIN 770 

79 


DIMENSION LSAVE(5C»203)» TLO ( 50 ) t NH OLO ( 15 ) 

MAIN 780 

ec 


DIMENSION SWPR0Ef(5O iTLDS (50 )#TLS(50) 

MAIN 790 

81 


DOUBLE PRECISION DP REC »OPR ECl # DPRE C 2 # DPREC 3 f OPRE C4 t 0PREC5 

MAIN SCO 

82 


DOUBLE PRECISION PIN T t Y E *P SN 5 f TO T t PRO B #PR 0 Bl- 

MAIN 810 

83 
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MAIN 320 
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MAIN 830 

85 

C 


MAIN 84G 

86 


Mi = 50 ^MAXIMUM NUMBER OF BLOCKS IN DIAGRAM* 

MAIN 850 

87 


M2 = 200 ^MAXIMUM NUMBER OF SUCCESS PATHS. 

MAIN 860 

88 


K3 r 14 ^MAXIMUM NUMBER OF INPUTS/OUTPUTS TO/FROM ONE BLOCK. 

MAIN 870 

89 


M4 = 15 ^MAXIMUM NUM3ER OF SENSE BLOCKS CONTROLLING STANDBY. 

MAIN 880 

SO 


M5 = 29 ^MAXIMUM NUMBER OF STANDBY BLOCKS CONTROLLED BY 1 SENSOR 

.MAIN 89C 

91 


Mo = 20 ^MAXIMUM NUMBER OF EQUIVALENT BLOCKS IN A SINGLE SET. 

MAIN 300 

32 


M7 = 20 a'MAXIMUM NUMBER OF EQUIVALENT BLOCK SETS. 

MAIN 910 

93 


Ll-Ml+i 

MAIN 920 

94 


L3-H3+1 

MAIN 930 

95 


L4i:M4 + Ml 

MAIN 940 
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L5rM5*U 

MAIN 950 
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LSr (2*M6)*1 

MAIN 960 
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LL6-M6-M 

MAIN 970 
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THE FOLLOWING SHOWS THE RELATIONSHIP OF THE ABOVE TO STORAGE. 

MAINIOCO 

1C2 
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CGMM0N/ALLSUB/Ml»K2#M3fM 4*M5f MS* M7 t LI t L3 #L4 tL5 • L6 

MAIN1020 

104 

C 

CCMMCN/LGSUE/IB(L3t2 *Mi ) t IIS (M 1 ) * I RE ( Ml t L5 )* ITEKP(L6 * M7 ) tNl t 

MAIN1C3C 

1C5 

c 

1 ISAVEIMI* M2 ) fISJC* JSUC *JSlf JSAVEC Ml* M2 ) * I JS 1 # T OT * PROS ( L 4 ) # 

MAIN104Q 

1C6 

c 

2 FRCtl(L4)fIf'RINTfNSTDBY(M4) tNSBYMX 

MAIN1G5C 

1C7 

c 

C0MM0N/SMSU8/NPUP ( M2) *N PUP MX # NS OUT (Ml) * NO OTMX • NBIN (Ml ) # N3INMX t 

MAIN106G 

ice 

c 

1 NBNUM( Ml ) *NtNM AX 

HAIN107C 

1C9 
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MAIN1080 

lie 

c 

NST - NSTEP 1 

MAIN109C 
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c 

DIMENSION RQ( Ml) *NP (Ml)* RT ( L4 *NST ) * PS NS ( NST ) * IS < Ml ) * T L C Ml ) 

MAINIICC 

112 

c 

DIMENSION IN (Ml ) tICUKMl) *JRB(L5) *NRB(L5) *KRe(L5) •KSAVECKl* M2) 

MAIMIID 

113 

c 

DIMENSION LSAVElMlf M2 > * IDR (Ml ) # TLO ( Ml ) »NHOLD ( L3 ) 

MAINX120 

114 

c 

DIMENSION SWPROBC Ml ) * T L DS ( Ml ) * TL S ( Ml ) 

MAIN113C 

115 
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MAIN1131 

116 

c 


MAIN1132 

117 

c 

THE variable NS13 SPECIFIES THE NUMBER OF *SIG NIFIC ANT ♦ FIGURES 

MAIN1133 

118 

c 

TO BE PRINTED FOR THE RELIABILITIES. NOTE THAT ♦SIGNIFICANT* 

MAIN1134 

119 

c 

FIGURES IS DEFINED AS THE NUJ^BER OF NON-NINES IN THE RELIABILITY 

MAIN1135 

12D 

c 

NUMBER. THUS* .33360* .385* AND. 56 ALL HAVE TWO ♦SIGNIFICANT* 

KAIN1136 

121 

c 

DIGITS. 

MAIN1137 

122 
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MI AN1138 

123 


N3IG-3 

MAIN1139 

124 

c 


MAIN114C 

125 

c 

THE ARRAY V HOLDS THE FORMAT FOR THE FINAL RELIABILITY PRINTING. 

MAIN1141 

126 

c 

THE ARRAY VK HOLDS THE ALPHA NUMBERS NEEDED TC CHANGE THE 

MAIN1142 

127 

c 

VARIABLE FORMAT V. 

MAIN1143 

128 

c 


MA3N1144 

129 


DIMENSION V(15)*VK(10) 

MAIN1145 

130 

c 


MAIN1146 

131 


DATA V/'(36HCR’**ELIA3I’ **LITY 0***F THE *f*SYSTEM’*» THRU ^* 

MAIN1147 

132 


i*TIME* f •FIC.E** *’3H HOU***RS = •** F* t* 10 *#*.** ’S* * M V 

MA2N1148 

133 


DATA VK/' 1» • •2**'3***4»**5***6«*'7***3***9’»*1Q*/ 

MAIN1149 

134 
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MAIN115C 

135 

c 

CLEAR VARIABLES TO ZERO AND SET CONSTANTS. 

MAIN116Q 

136 
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MAIN117C 

137 


IPRINT^G 

MAIN1175 

138 


393 LASTED 

MA1N118Q 

139 


NOPRNT=i 

MAIN1190 

140 


DC IS I=1*L3 

MAIN120G 

141 


00 19 J=l*2 

MAINX210 

142 


DO 19 K=1*M1 

MAIN1220 

143 


19 I3( I* J*K )-D 

MAIN123Q 

144 


DO 6 I=l*Mi 

MAIN124C 

145 


IS(I)--C 

MAIN125C 

146 


IIS(I)=D 

MAIN126C 

147 


DO 6 K=i*L5 

MAIN1270 

148 


6 IRB(I*K):D 

MAIN128G 

143 
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MAIN1290 

150 


397 REWIND 1C 

MAIN13CC 

151 


NS3YMX=Q 

MAIN132C 
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152 


KTrC 

MAIN133C 

153 


.\6H0L0=a 

MAIN1340 

154 


NSTEF=1CG 

MAIN135C 

155 


JJJ=0 

MAIN136G 

156 


KRBl=i 

MAIN1370 

157 


NR31-Q 

MAIN138C 

158 


CO 1C I-ltL4 

MAIN139G 

159 


PROSMUC.OO^C 

MAIN14C0 

16C 


PR0B1(I)-C*CL+D 

MAIN1410 

161 


00 10 L-l»101 

MAIN1420 

162 


RTd fLl-0. 

MAIN143C 

183 


1C PSNS( D-O.OO+G 

MAIN1440 

164 


uC 15 I-1»M4 

MAIN145G 

165 


15 NSTDaY(I)=Q 

MAIN1460 

166 


DO 16 I-ltL3 

MAIN147C 

167 


IN( I)=0 

MAIN148Q 

168 


16 ICUT(I)=C 

MAIN1490 

169 


OC 20 I-1»M1 

MAIN1500 

170 


SWPRCBd J^l.C 

MAIN151C 

171 


TLOS( I)=u. 

MAIN1520 

172 


TLS( I)=0. 

MAIN1530 

173 


TL( I)=0. 

MAIN1540 

174 


TLD(I)=0. 

MAIN155G 

175 


Ru( I)=C. 

MAIN156Q 

176 


NPdI-C 

MAIN157G 

177 


OC 17 K-i»L5 

MAIN1580 

178 


NRB( K)=G 

MAIN159C 

179 


17 KR3(K)=Q 

KAINIGOO 

18C 


OC 2C K=1»M2 

MAIN161G 

181 


JSAVEt If K)-0 

MAIN1620 

162 


2G ISAVC(IfK)=0 

MAIN163G 

183 

C 


HAIN1640 

184 


IF<LAST.LG.7 ) GO TO 30CC 

MAIN165C 

185 

C 


MAIN1660 

166 

C 

READ INPUTS AND OUTPUTS TC BLOCK NUMBERS. 

MA1N1670 

187 

c 


MAIN1680 

168 


iC99 READ (5fl2CG »ERR::3CC#ENC-3CC ) NB t IN » IC UT t LAS T 

MAIN1&9C 

189 


12QC FORMAK 3iI2fi7Xf ID 

MAIN1700 

190 

c 


MAIN171C 

131 

c 

PRINT DATA CAROS. 

MAIN1720 

192 

c 


MAIN173C 

193 


IF( IPRINT,N£.4) go TO 1210 

MAIN1740 

1S4 


PRINT 12C2# NBflfJtlOUTtLAST 

MAIN175C 

195 


1202 FORMAK* DATA CAR0«31I3f» LAST-’IZ) 

MAIN1760 

1S6 

c 


MAIN177C 

197 

c 

NBHOLD - TEMPORARY HOLD FOR THE HIGHEST NUMBER BLOC K ( N BMAX ) • 

MAIN178Q 

198 

c 


MAIN179C 

193 


121C IF( NBHOLD-NB) 1C2G f 1021 f 1 022 

MAINldOO 

2C0 


1G2G NBHOLDrNB 

MAIN1810 

201 


GO TO 1022 

MAIN1820 

2C2 


IC'21 PRINT lC23f Nb 

MAIN183G 

203 


1023 FORMATt dERROR IN MAIN PROGRAM AT 102 3 . V * OMORE THAN ONE BLOCK NUMMAIN1840 

2C4 


IbER -U4) 

MAIN185C 

205 


GO TO 1C68 

MAIN1860 

206 

c 


MAIN1861 

2C7 

c 

REREAD CARO IF IPRINT IS 3EING SPECIFIED 

MAIN1862 

2C8 

c 


MAIN1863 

209 


30C READ(Cf 301fERR::3Q2f EN0 = 3G2) I#J 

MAIN1364 

21G 


3C1 FCRMAT(7XiIl tGXtID 

MAIN1865 

211 


IF( (I.6T.4).CR.(J.GT.3n GO TO 3G2 

MAIN1866 

212 


IPRINT=I 

MAIN1867 

213 

214 

215 


IF<J. EO.C) j=3 
NSIG = J 

PRINT 3G3#IPRINT#NSIG 

MAIN1368 

216 


303 FORMAT! *GIPRINT OPTION SET AT»I4** NSI6 OPTION SET AT»I4/1H1I 

MAIN1869 

217 


GO TO 1039 

MAIN1870 

218 


3C2 PRINT 304 

MAIN1871 

219 


3G4 FORMAT! •GERRCR** ATTEMPT TO SPECIFY IPRINT OR NSIGt BUT SPECIFIED 

MAIN1872 

220 


ICUT OF RANGE.*/* CARD IGNORED*) 

MAIN1673 

221 


GO TO iC33 

MAIN1874 

222 

c 


MAIN1875 

223 

c 

GENERATE ARRAY IB WITH IBdflfNB) = NUMBER OF INPUTS TO NBt 

MAlNld80 

224 

c 

IBCltZfNB) = NUMBER OF OUTPUTS FROM NE» 

MAIN189C 

225 

c 

THE REST OF THE 13 ARRAY CONSISTING OF INPUTS AND OUTPUTS TO NB. 

MAIN190Q 

226 

c 


MAIN191C 

227 


1022 DC IClu I-lfL3 

MAIN1920 
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226 


IFdNdJ.OG.U GO TO 1C il 

MAIN1S30 


229 

lOlG 

CCNTIKUE 

MAIN1S40 


23C 

IGli 

INl-1-1 

MAINX950 


231 


DC 1G12 I-lfL3 

MAIN1960 


232 


IF(IGUT(I)«EG*0) GO TO 1D13 

MAIN1970 


233 

1C12 

CONTINUE 

MAIN1980 


234 

1C13 

ICUT1=I-1 

MAIN199C 


235 


IB( Idf N3) = lNi 

MAIN20CO 


236 


IbatZtNB) = lOUTl 

MAIN201D 


237 


00 1C14 I=IfINl 

MAIN2020 


238 

1014 

Ib(I+ltl#NE)=IN(I) 

MAIN2C3C 


239 


DC 1C15 I^lflCUTl 

MAIN2C4D 


240 

1C15 

lL(l4l,2fNB)::I0UT< I) 

MAIN2C5C 


241 


IF(LAST.NE«9> GO TO 1099 

MAIN2060 


242 

C 


MAIN2C7C 


243 

c 

CHECK SUCCESS BLOCK OUTPUT - ZERO. 

MAIN2C80 


244 

c 


MAIN209G 


245 


DC 1C25 I-lfL3 

MAXN21C0 


246 


IF (ICUTd J.EO.C ) GC TO 1G25 

MAIN2110 


247 


PRINT lG24f lOUT(I) 

MAIN2120 


248 

1C24 

FORMAT( *1ERRCR IN MAIN PROGRAM AT 1C24 . • / «CS LCCESS BLOCK HAS OUT PUMA IN213 0 


249 


LT -d4) 

MAIN2140 


250 


GO TO 1C66 

MAIN2150 


251 

1C25 

CONTINUE 

MAIN2160 


252 

c 


MAIN217C 


253 

c 

I5UC = SUCCESS SLOCK NU«8ER(BASE OF PROBABILITY TREE TO BE GENER- 

MAIN2180 


254 

c 

ATED BY CALLING SUBROUTINE TREE). 

MAIN219C 


255 

c 


NAIN2200 


256 

1017 

ISUC=NB 

MAIN2210 


257 

C 


MAIN2220 


258 

c 

J5UC AND IPRINT ARE USED IN SUBROUTINE TREE TO CONTROL PRINTING. 

MAIN2230 


259 

c 


MAIN2240 


260 


JSUC-NB 

MAIN2250 


261 

c 


MAIN226C 


262 

c 

NBMAX r The HIGHEST BLOCK NUMBER OF BLOCKS 1-50/1-Ml* 

MAIN2270 


2£3 

c 


MAIN228C 


264 


NEMAX^NBHOLO 

MAIN229C 


265 

c 


MAIN230G 


266 

c 

PRINT RELIABILITY BLOCK DIAGRAM. 

MAIN2310 


267 

c 


MAIN232C 


263 


PRINT 632 

MA1N233C 


269 

632 

FORMAT! ^RELIABILITY BLOCK DIAGRAM*) 

MAIN234C 


270 


DO IS I-ifMl 

MAIN235C 


271 


IF(IBd»ItI ) )21#22f21 

MAIN236C 


272 

22 

lF(iad»2»I))21»18.21 

MAIN2370 


273 

21 

KdE(l*l*I) + l 

MAIN238C 


274 


IF(K.LE.I) GC TO 4C 

MAIN2381 


275 


PRINT 202*1 • (IE(Jfl»I)»J-2fK » 

MAIN2382 


276 

2C2 

FORMAT! lHG5HaL0CKl3f 3X*5 HINPUT14 ! 13 *iX ) ) 

MAIN2383 


277 


GC TC 41 

MAIN2384 


278 

4C 

PRINT 2G5*I 

M AIN2385 


279 

205 

FORM AT! 1HC5HEL0CKI3 *5X* 5HINP UT ) 

MAIN2386 


280 

41 


MAIN2387 


261 


IF!K.LE.l) GC TO 42 

MAIN2388 


232 


PRINT 2G3*!Ib!J*2#I)#J=:2*K) 

MAIN2383 


283 

2C3 

FORM AT! 11X*6HCUTPUT *14! I3*1X )) 

MAIN239C 


234 


GC TO 13 

MAIN2391 


285 

42 

PRINT 2G6 

MAIN2392 


286 

2C6 

FORMAT! 1 IX* 6H OUTPUT) 

MAIN2393 


267 

18 

CONTINUE 

MAIN24DC 


288 

C 


MAIN241G 


289 

C 

READ SENSE BLOCK AND STANDBY BLOCKS OF THAT SENSE BLOCK. 

MAIN242C 


290 

c 


MAIN2430 


291 

c 

IIS ARRAY holds SENSE BLOCKS OF ORIGINAL TREE. 

MAIN244C 


292 

c 

IS ARRAY HOLDS SENSE BLOCKS OF ORIGINAL TREE AND VARIES. 

MAIN2450 


293 

c 

NNSR = NUMBER OF SENSE BLOCKS IN ORIGINAL TREE. 

MAIN246C 


294 

c 

NSR ^ NUMBER OF SENSE BLOCKS. IN ORIGINAL TREE AND VARIES. 

MAIN247G 


2S5 

c 

IRB ARRAY HOLDS STANDBY BLOCKS CONTROLLED BY THEIR SENSE BLOCKS. 

MAIN248C 


296 

c 

IRB(NUM5£R OF SENSE BLOCK* D - MAXIMUM NUMBER OF STANDBY BLOCKS 

MAIN249C 


297 

c 

OF THAT SENSE BLOCK. 

MAIN250C 


293 

c 

IR8!NUMSER OF SENSE BLOCK* No) - STANDBY BLOCK NUMBER CONTROLLED 

MAIN2510 


299 

c 

BY THAT SENSE BLOCK. 

MAIN252C 


3Q0 

c 


MAiN253C 


3C1 


READ 12G0t NSR 

MAIN254C 


302 


NNSR- NSR 

MAIN2550 


303 


IF!NSR.EG.Q) GO TO 84C 

MAIN256C 
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3C4 


DO 2CC1 I=lfNSR 

MAIN2570 

3Q5 


READ 12CG N5*JRB 

MAIN258C 

306 


irs(i)=Ns 

MAIN259Q 

3G7 


IS(I)-NS 

MAIN26CC 

308 


DO 2002 J=1»M5 

MAIN2610 

309 


IREC NS» J + 1) = JRB( J) 

MAIN262C 

310 


IF ( jR3t J) .Ea.C) GO TO 2003 

MAIN263C 

311 

2CC2 

CONTINUE 

MA3N2640 

312 

2003 

IR8(NSf i )>U-1 

MAIN2650 

313 

2C01 

CONTINUE 

MAIN266C 

314 


DO 2020 I-1»NSR 

MAIN2670 

315 


INS=IS(I) 

MAIN268C 

316 

C 


MAIN2690 

317 

c 

STANDBY 3L0CKS(NH0LDJ CONTROLLED BY A SENSE cLOCKIINS) MUST 

F0LLCWMAIN270C 

318 

c 

ALL OTHER INPUT BLOCK NUMBERS IN IB ( N8 » 1 1 1 NS J • 

MAIN2710 

319 

c 


MAIN272C 

320 


DC 2200 J=1*L3 

MAIN273Q 

321 


NHOLD* J)::C 

MAIN274C 

322 

2200 

CONTINUE 

MAIN275C 

323 


NUM=C 

MAIN276C 

324 


KEND-IRBt ZNSf l) + i 

MAIN277Q 

325 

c 


MAIN278C 

326 


DO 2250 L=2fL3 

MAIN279Q 

327 


IFaE(Lti#INS).E3.C) 60 TC 2260 

MAIN2800 

328 


DO 2240 K=2fK£ND 

MAIN2810 

329 


IF{IB(LtltINS).NE.lRB(INSfK) > GO TO 2240 

MAIN262C 

330 


NUMrNUM+1 

MAIN283Q 

331 


NHOLD(NUM J-IP.Bt INSiK ) 

MAIN284C 

332 


I6(L#1»INS)=C 

MAIN285C 

333 


60 TC 2250 

MAIN2860 

334 

2240 

CONTINUE 

MAIN287Q 

335 

2250 

CONTINUE 

MAIN288C 

336 

2260 

IF( NUM.NE.O) GO TO 2300 

MAIN2890 

337 


PRINT 227C» INStNUMfINS 

MAIN290C 

338 

2270 

FORMAT! *1ERR0R IN MAIN PROGRAM AT 2270 . V * CERROR DURING SENSE BL0CMAIN291C 

339 

iK»I3.’ INPUT REARRANGE TO STCRO ACTIVE BLOCKS BEFORE STANDBY 

BL0CKMAIN292C 

340 

2S IN I3( ’I3»* »lf • I3f M .• ) 

MAIN2930 

341 


GO TC 1066 

MAIN294C 

342 

C 


MAIN29SC 

343 

C 

SHIFT IB(NE»iflNS> TC LEFT ELIMINATING ZEROES. 

MAIN296C 

344 

C 


MAIN2970 

345 

23GC 

ILEFT=G 

MAIN238G 

346 


00 2350 J=2tL3 

MAIN2990 

347 


IF(Ib( J»lfINS).NE.C> 60 TO 2330 

MAIN30DG 

343 


IL£FT=ILEFT*1 

M AIN3010 

349 


GO TC 235D 

MA3N3C2G 

350 

2330 

IF( ILEFT.EQ.G ) 30 TO 2350 

M AIN303Q 

351 

C 


MAIN3040 

352 

C 

MOVE N3 TO LEFT. 

MAIN3G5C 

353 

C 


MAIN3C6C 

354 


J--ILEFT 

MAIN307Q 

355 


IE(M#l»INS)=IB(Jtl»lNS) 

MAIN308G 

356 

2350 

CONTINUE 

M AIN3C9Q 

357 

C 


MAIN31CC 

358 

C 

ADO N3 IN NHOLD AT END. 

MAIN3110 

559 

c 


MAIN312C 

330 


KoEG=L3-Il£FT+1 

MAIN3130 

361 


M = C 

MAIN314C 

362 


DC 2370 K-K8EG*L3 

MAIN3150 

363 


M::M+I 

MAIN316C 

364 


I3(K#1»INS)=NH0L0!M> 

MAIN317Q 

365 

2370 

CONTINUE 

MAIN318C 

366 

C 


MAIN3190 

367 

c 

PRINT SENSE BLOCK AND STANDBY BLOCKS OF THAT SENSE BLOCK. 

MAIN320C 

368 

c 


MAIN321C 

369 


DO 2025 K=2tL5 

MAIN3211 

370 


IF( IRB( INS fK) .EQ.O) GO TO 2026 

MAIN3212 

371 

2025 

CONTINUE 

MAIN3213 

372 

2026 

KrK-1 

MAIN3214 

373 


PRINT 2G21f INSf (IRB(lNSiJ)i J=2tK> 

MAIN322C 

374 

2G21 

FORMAT! •CSENSE SWITCH*I3»* CONTROLS THE FOLLOWING SLOCKS* # 2 !/lH t MAIN3230 

375 


15CI4 ) ) 

MAIN3231 

376 

2020 

CONTINUE 

MAIN3239 

377 

c 


MAIN324C 

378 

c 

READ AND GENERATE ITEMP ARRAY CF EQUIVALENT SLOCKS. 

MAIN3250 

379 

c 

z MAXIMUM NUMBER OF EQUIVALENT BLOCK SETS. 

KAIN326C 


6 


JPL Technical Memorandum 33-513 



38C 

c 


MAIN3270 

381 

640 

RCAD 12DD» Ni 

MAIN3280 

332 


IFCNI. 6:3.0) GO TO 2040 

MAIN329G 

383 

141 

DO 12.5 I = 1*N1 

MAIN33CC 

384 

125 

READ l24#(ITEMP(>J»I)f J = 2#LL6) 

MAIN3310 

385 

124 

FORM AT(4CI2 ) 

M A IN 332 C 

386 

C 


MAIN333G 

3£7 

C 

PRINT equivalent BLOCKS. 

MAIN334C 

388 

C 

MAIN335Q 

389 

140 

PRINT 905 

MAIN336C 

390 

9G5 

FORMAT( lH017HEQUiVALENT BLOCKS) 

M AIN3370 

391 

DC 123 I^lfNl 

MAIN338C 

392 


PRINT 3C6f (ITEMP( J#D# J-2#LL6> 

MAIN339C 

393 

906 

FORMAT ( •C*4DI3 ) 

MAIN340C 

394 


CG 123 K=2fLL5 

MAIN3410 

395 

123 

ITEMF ( K + M6 tl J 1 TEMP t K • I) 

MAIN342C 

396 

2C4C 

CONTINUE 

MAIN3430 

397 

C 


MAIN344C 

398 

c 

RELIABILITY BLOCK DIAGRAM INPUT CHECKED. 

MAIN3450 

399 

400 

c 

c 

MAIN346C 

MAIN3470 

401 

c 

GENERATE CRItINAL TREE. 

MAIN348C 

4C2 

c 

ISAVE HOLDS THE PROBABILITY TREE SUCCESS PATHS 

MAIN349Q 

4C3 

c 

JSI = NUMBER OF SUCCESS PATHS. 

MAIN3500 

4G4 

c 

MAIN351Q 

4C5 

30GD 

CALL TREE($99£) 

MAIN352C 

406 

c 

MAIN3530 

407 

C 

DUPLICATE ORIGINAL TREE INTO JSAVE. 

MAIN354D 

4C8 

c 

MAIN3550 

4C9 


CALL DUPTRE(JSAVEflJSli ISAVE t JSI) 

MA1N356C 

41C 

c 


MAIN357Q 

411 

c 

IDIST DETERMINES WHAT FAILURE DISTRIBUTION IS USED < 1 ::EXPONEN TIAL ) 

MAIN358C 

412 

c 

TTOT - MISSION TIME. IPRINT 1 = PRINT ALL TREESf 

MAIN3590 

413 

c 

0 = SUPPRESS PRINTING OF ALL TREES EXCEPT THE ORIGINAL TREE. 

MAIN360C 

414 

c 

MAIN361Q 

415 


READ 40CCt IDlSTt TTOT 

MAIN362G 

418 

4CCC 

FORMAT! I2f £12 .7) 

MAIN363C 

417 

c 


KAlN3e4G 

418 

c 

READ PARAMETERS OF FAILURE DISTRIBUTION. 

MAIN3650 

419 

c 

RC IS SET = l.C IF READ IN AS BLANK OR ZERO. 

HAIN3660 

42C 

c 

MAIN3670 

421 


GO TO ( 4100 » 42 CO »4 300 t4 t4 SCO ) t IDIST 

MAIN368C 

422 

c 


MAIN3690 

423 

c 

READ ACTIVE PARAMETERS. 

MAIN37DD 

424 

c 


MAIN3710 

425 

41CC 

READ 4 1C If NBiTLAMBDf RLltLAST 

MAIN372C 

426 

4101 

FORMAT! I2t£12 .7f Fl0.7»55Xf ID 

MAIN373Q 

427 


IF!RDl)41C9f4110t41C3 

MAIN374C 

428 

41G9 

TL! NB )=TLAMBD 

MAIN3750 

429 


RC!NE)=R03 

MAIN376C 

430 


GC TO 412C 

MAIN3770 

431 

4130 

R D ! N B ) = 1 . 0 

MAIN3780 

432 


TL! NB )=TLAM3D 

MAIN379Q 

433 

4120 

IF (LAST.LT.5 > 60 TC 4 1C C 

HAIN3800 

434 

C 


HAIN381C 

455 

C 

READ DORMANCY PARAMETERS. 

MAIN3820 

436 

C 

DFACT : DORMANCY FACTOR TO BE MULTIPLIED BY STANDBY ACTIVE LAMBOASM AIN383Q 

437 

C 

IF IDUMMY NON ZERO AND DFACT = 0 - READ LAMBDA DORMANT = TLDCNBl. 

HAIN384C 

433 

c 


HAIN3850 

439 


IF(NSR.E3.G) GC TO 471D 

MAIN386C 

440 


READ 4101f lOJKMYf DFACT 

MAIN3870 

441 


IF ilDUMMY.EG.C ) GC TO 416D 

HAIN3880 

442 


IF!0FACT.GT.D.) go to 416C 

MAIN3890 

443 

4150 

READ 4iClf NBf TLAMBD-t DUMMY? LAST 

MAIN39C0 

444 

TLD !Nt » = TL AMBD 

MAIN391Q 

445 


IF !LAST.LT.6 ) GC TO 415 C 

MA2N392C 

446 


SO TO 4190 

MAIN3930 

447 

4160 

DC 418C J-lfNSR 

MAIN394C 

44 3 


KSrIS! J) 

MAIN3950 

449 


KEND::IRb! NSf D ♦! 

MAIN396C 

450 


DC 4170 K:::2fKEND 

MAIN3970 

451 


NL-IRB ( NS fK ) 

MAIN398C 

452 


TLD!N5)=DFACT*TL!Ne) 

M AIN3990 

453 

4170 

CONTINUE 

MAIN4CCC 

454 

4180 

CONTINUE 

MAIN4D10 

455 

C 


HAIN4C2C 
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45S 

C 



MAIN4030 

457 

C 

READ SWITCHING OPTIONS AND DATA FOR EACH SENSE BLOCK. 

MAIN404C 

458 

C 



MAIN4050 

459 

c 

0 PERFECT SWITCHING( PROBABILITY SWITCH WORKS 

EQUALS l.Cl 

MAIN4060 

460 

c 

1 - CONSTANT PROBABILITY THAT SWITCH WORKS. 


MAIN4070 

461 

c 

2 = CORKANT FAILURE RATE FOR SWITCH. 


MAIN408C 

462 

c 

3 “ DORMANT AND ACTIVE FAILURE RATE FOR SWITCH. 


MAIN409C 

463 

c 

4 = NOT available 


MA3N4100 

454 

c 

5 - NOT AVAILABLE 


MAIN4110 

465 

c 



MAIN412C 

466 

4190 

DO 4700 J-1*NSR 


MAIN4130 

467 


READ 41Cif NEf TLAMBOtSPROBfISCPT 


MAIN4140 

468 


IF( ISOPT.EQ.C ) GO TO 4600 


MAIN415Q 

469 


60 TO <4&lGt462Ct463Ct4£4Gt465C)* ISCPT 


MA1N416C 

470 

46C0 

SWPR03( NB )-l. C 


MAIN4170 

471 


TLDS (NB)rO.C 


MA2N418C 

472 


TLS(N6)::G.O 


MAIN4190 

473 


60 TO 4680 


MAIN420C 

474 

4610 

SWPR03( NS) ^SPROB 


MAIN4210 

475 


TEDS (NB)=D.C 


MAIN422C 

476 


TLS(Nb)-O.C 


MAIN4230 

477 


60 TO 4660 


MAIN4240 

473 

4620 

SWPR03(Ne)=l.G 


MAIN4250 

479 


TLDS(NB)=TLAMBD 


MAIN426C 

480 


TLStNB)::Q.C 


MAIN427Q 

461 


6C T.C 4680 


MAIN428C 

432 

4630 

SWPR03(N5)=1.G 


MAIN4290 

463 


TLDS(NB)=TLAMED 


MAIN4300 

434 


READ 41U1» NBfTLAMBD 


MAIN4310 

465 


TLS( Nb)=TLAMBD 


MAIN4320 

486 


GO TO 4630 


MAIN4330 

467 

484C 

C0NTI^4UE 


MAIN434D 

438 


GO TO 4530 


MAIN4350 

489 

465C 

CONTINUE 


MAIN436C 

490 

4680 

TLD(N& )=G.G 


MAIN4370 

491 


TL(NE) -O.C 


MAIN438C 

492 


R01NB)=1.Q 


MAIN4390 

493 

47CC 

CONTINUE 


HAIN4400 

494 

C 



MAIN4410 

495 

C 

READ LAST(8C)» LAST = 7* 8t CR 9. 


MAIN442C 

496 

C 

7 - RECALCULATE WITH NEW PARAMETERS^ 8 : NEW OIAGRAMt 9 - END. 

MAIN4430 

497 

C 



HAIN4440 

498 

4710 

READ 27C7. LAST 


MAIN445Q 

499 


GO TO 4600 


MAIN446C 

SCO 

C 



MAIN4470 

501 

42DC 

CONTINUE 


MAIN448C 

5C2 

430C 

CONTINUE 


MAIN4490 

503 

44CC 

CONTINUE 


MAIN45GC 

504 

45CC 

CONTINUE 


MAIN451G 

505 

C 



MAIN4520 

506 

C 

42CC»43Q0f 44C0»450C AVAILABLE IF YOU WANT TO USE 

OTHER THAN 

MAIN453C 

5C7 

C 

EXPONENTIAL US T R IBU T 10 N . ADD APPROPRIATE STATMENTS FOR NEW 

MAIN454C 

508 

C 

DISTRIBUTION. 


MAIN4550 

509 

C 



MAIN4560 

510 

4300 

IF(NSR.EQ.a) GO TO 8CCC 


MAIN4570 

511 

C 



MAIN458C 

512 

c 

FIND A SENSE TREE WITH NC OTHER SENSE BLOCKS IN 

IT. 

MAIN4590 

513 

c 



MAIN460C 

514 

7CCC 

CONTINUE 


MAIN4610 

515 


DO 7010 JJ=1#M1 


MAIN462C 

516 

7010 

NP( JJ)=C 


MAIN463C 

617 


INDEPro 


MAIN464Q 

519 


M-1 


MAIN4650 

519 


J-1 


MA1N466C 

520 


I-l 


MAIN4670 

521 

7011 

IF(JSAVE(If J).N£.0) GO TC 7CC1 


MAIN4680 

5 22 


IF( I-l) 7002# 7003 »7CG2 


MAIN4690 

523 

7001 

Krl 


MAIN470G 

524 

7004 



MAIN4710 

525 

7005 

CONTINUE 


MAIN472C 

526 


IF(NP<L).EQ.IS(K> ) GO TO 7005 


M AIN473C 

527 


IF (L.EQ.M) 60 TO 7CC7 


KAIN474C 

528 


L=L+1 


M AIN4750 

529 


GO TO 7005 


MAIN4760 

530 

7CG7 

IF( A8S< JSAVE( I» J) ) .E3.IS (K ) ) GO TO 7008 


MAIN477C 
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531 

7006 

532 

533 

534 

7008 

535 

536 

537 

538 

539 

540 

7009 

541 

542 

543 

7002 

544 

545 
54S 
547 

7003 

548 

549 

550 

20QQ1 

551 

552 

553 

20C02 

554 

C 

555 

C 

556 

C 

557 

c 

558 

c 

559 

c 

560 

c 

561 

c 

562 

c 

563 

c 

564 

7030 

565 

566 

567 

568 
£69 
570 

7C5G 

571 

C 

572 

573 

C 

574 

c 

575 

c 

576 

577 

578 

579 

580 

c 

581 

c 

532 

c 

583 

c 

584 

c 

585 

c 

586 

c 

587 

7200 

588 

c 

569 

590 

c 

591 

c 

592 

c 

593 

594 

595 

c 

596 

597 

c 

538 

c 

599 

c 

600 

601 

602 

603 

604 

605 

7400 


IF(K.tQ.NSR) GO TO 7DG9 
K-K+1 

GC TO 7CCH 

IF( ISC K)*£G.INO£P) 6C TO 7009 

NPC M-1 )-INOEl? 

INDEF-IS(K) 

IF( IPRINT.NE.4) GO TO 7003 
PRINT 2CCC1# INDEP 
IF( I.EQ.Ml ) GC TC 7CG2 
1 = 1+1 

GO TO 7011 

IFCJ.EG. M2) GC TC 7D03 

1=1 

J=J+1 

GO TO 7011 
ISUC=INDEF 

IFC IPRINT.NE.4) GO TO 7030 
PRINT 200C1* ISUC 
FORMAT CI3) 

PRINT 2CDC2f ( I S ( I ) 1 1=1 .Ml ) 

PRINT 2GCC2. ( NP ( I ) t 1= 1. Ml ) 
F0RMATC1HC2 5 CI3*1X)/1H 25 (13 . IX)) 


MAIN478C 

MAIN4790 

MAIN48CC 

MAIN4810 

MAIN482C 

MAIN483Q 

MA1N484G 

MAIN4850 

HAIN486C 

MAIN4870 

MAIN488C 

MAIN4890 

MAIN49DC 

MAIN4910 

MAIN492C 

MAIN493Q 

MAIN494C 

MAIN4950 

MA1N4SSC 

MAIN4970 

MAIN498C 

MAIN4990 

MAIN500C 

MAIN501G 


DETERMINE WHAT TYPE OF STANDBY IS INVOLVED. 


HAIN502G 
MAIN5030 
MAIN5C4C 
MAIN505Q 

IF SENSE BLOCK CONTROLS ONE STANDBY BLOCK. SET ISTDBY = It SEE IF MAIN5D6C 
THAT STANDBY BLOCK HAS AN INPUT. IF NO INPUT GO TO 7700. IF IT MAIN5070 
HAS AN INPUT GC TO 75CD . IF SENSE BLOCK CONTROLS MORE THAN ONE MAIN508C 
STANDBY BLOCK. SET IST0BY=2 AND CALL SUBROUTINE STDBY2. MAIN5Q90 

MAIN51DC 

IF( IR3( INDEP. D.GT.l) GO TO 7C5C MAIN511C 
ISTDBY=1 MAIN512C 
IRBB=IR3( IN0LP.2) MAIN5130 
ITEST=IE(l.lfIRBB) MAIN514C 
1F( ITEST.EQ.Q I GO TO 77D0 MAIN5150 
GC TC 75GC MAIN516C 
IST0BY=2 MAIN5170 


CALL ST0BY2U»393. INDEP ) 

DUPLICATE STANDBY TREE RETURNED IN ISAVE INTO LSAVE. 

CALL OUPTREILSAVEf LJSlfl SAVE. JSl) 

IF(IFRINT.LT.3) GO TO 720C 
CALL TRPRNT(2.JSAVE.IJS1 f 1 ) 

CALL TRPRNTC 4.LSAVE.LJS1. 2 ) 

GENERATE FULL STANDBY TREE. 

SUBROUTINE TRENUM STORES IN NENUM ALL THE DIFFERENT ABSOLUTE 
BLOCK NUMBERS IN THE STANDBY TREE. NUMOUT ZEROES OUT THESE 
BLOCK NUMBERS. 

CALL TRENUM(LSAVE.LJSl) 

ISUC=IR6(INDEP.2 ) 

IF STANDBY BLOCK HAS NO INPUT. 60 TO 7500. 

ITEST=IB(1.1.ISUC) 

IF( ITEST.EQ.C ) GO TO 7500 

CALL TREE( $993) 

ZERO OUT STANDBY BLOCKS IN STANDBY TREE. 

CALL NUM0UT($393. ISAVE. JSl) 

CALL SFLEFT ( $998 .ISAVE. JSl) 

CALL 0UP0'JT($9S3. ISAVE. JSl) 

IF(IFRINT.LT.3) 60 TO 74C0 
CALL TRPRNTd .ISAVE. JSl. 3 ) 

CALL TRENUMC ISAVE. JSl) 


MAIN518C 

MAIN5190 

MAIN52DC 

MAIN5210 

MAIN522C 

MAIN523C 

MAIN5240 

MAIN5250 

HA1N526C 

MAIN5270 

MAIN528C 

MAIN5290 

VALUEMAIN530C 

MAIN5310 

MAIN5320 

MAIN533C 

MAIN534Q 

MAIN535Q 

MAIN536C 

MAIN537Q 

MAIN538C 

MAIN539Q 

HAIN540C 

MAIN5410 

MAIN542C 

MAIN5430 

MAIN544C 

MAIN5450 

MAIN546C 

MAINS470 

MAIN548C 

MAIN5490 

MAIN55CC 

MAIN5510 

MAIN552C 
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&C6 

C 






HAIN553C 

6C7 

C 






MAIN554C 

608 

C CALCULATE THE PROSABIL IT Y (P I TREE) THAT YOU 

NEED 

THE 

STANDBY 

TREE 

.MAIN555C 

609 

C 






MAIN556C 

610 

C 






MAIN5570 

611 

7500 ISUC=INOEP 






MAIN558G 

512 

CALL TREE! $993) 






MAIN55SG 

613 

C 






MAIN56CC 

614 

C DUPLICATE SENSE TREE INTO KSAVE. 






MAIN561Q 

615 

C 






MAIN562G 

616 

CALL DUPTREtK SAVE fKJSlfl SAVE# JSl) 






MAIN5630 

617 

C 






MAIN564C 

618 

C SUBTRACT FULL STANDBY TREE FROM SENSE 

TREE 

GIVING 

PI 

TREE. 


MAIN5650 

613 

C 






MAIN5S6C 

620 

GO TO (7515»7513) .ISTOBY 






MAIN5670 

621 

C 






MAIN568C 

622 

C GENERATE FULL STANDBY TREE. 






MAIN569C 

623 

C 






MAIN57CC 

624 

7515 ISUC=IR3( INDLPf 2) 






MAIN5710 

625 

CALL TREE($996) 






MAIN572C 

626 

C 






MAIN573Q 

627 

CALL TRENUMdSAVEt USD 






MAIN574C 

523 

C 






MAIN5750 

629 

C INITIALIZE 






MAIN5760 

630 

C 






MAIN577G 

631 

7513 DC 7520 I=1»M1 






MAIN578C 

632 

NdOUK D-C 






MAIN5790 

633 

7520 CONTINUE 






MAIN58GC 

634 

C 






MAIN5810 

635 

C ELIMINATE PATHS IN SENSE TREE IF THEY 

CONTAIN THE 

STANDBY BLOCKS. 

MAIN582C 

636 

C 






MAIN5830 

637 

NEOTMXriReCINDEPfl) 






MAIN584C 

638 

DO 7530 J-l»NaOTHX 






MAIN585Q 

639 

NBOUK J)-IRB (INDEP# J^l) 






MAIN586C 

640 

7530 CONTINUE 






MAIN5370 

641 

CALL PATHCT($998fKSAVE»KJSl) 






MAIN588C 

642 

IF( IPRINT.LT. 3) GO TO 7535 






MAIN5890 

843 

CALL TRPRNT(3fKSAVEf KJSl* 4 ) 






MAIN59CC 

644 

C 






MAIN5910 

645 

7535 NENM AX=NBNMAX+1 






MAIN5S2C 

646 

NSNUM( N3NMAX) -INOEP 






MAIN5930 

647 

CALL NUMCUT($S98tKSAVEtKJSl) 






MAIN594G 

648 

CALL SFLEFTtiSga# KSAVEfK JSl) 






MAIN5950 

649 

CALL DUPCUT($SS8fKSAVEtKJSl) 






MAIN596C 

650 

IF( IPRINT.LT. 3) GO TO 7538 






MAIN5970 

651 

CALL TRPRNTI 3tKSAVE#KJSl# 5 ) 






MAIN5S8C 

652 

C 






MAIN599C 

653 

C CHECK PI TREE FOR STANDBY TREES REPLACED BY 

STANDBY 

BL0CK(51 

“65 )/ 

MAIN6C0Q 

654 

C <L1-L4). 






MAIN6010 

655 

C 






MAIN602C 

656 

7538 IF (NSBYMX.LE.l) GO TO 7545 






MAIN6030 

657 

REWIND ID 






MAIN6C40 

658 

DO 7540 N=i#NSBYMX 






MAINS050 

659 

READ(IC) NSSfJSl 






MAIN&C6C 

680 

DO 7537 I-1#M2 






MAIN607C 

661 

RE AD (10) (ISAVE( J»I) #J-ifMl) 






MAIN6C8G 

662 

7537 CONTINUE 






MAIN609C 

663 

CALL TRINTR(S99 8 t KSAVE fKJ Sit ISA VE» JSl 

tNSTDBYtNSBYMX* 

NSS) 


MA1N61CC 

6 64 

7540 CONTINUE 






MAIN6110 

665 

C 






MAIN612C 

866 

C DUPLICATE KSAVE INTO ISAVE. 






MAIN6130 

667 

C 






HAIN614C 

663 

7545 CALL DUPTRE ( I 5 AVE • JS 1 #K SA VEt K JS 1 ) 






MAIN6150 

669 

C 






MAIN616C 

67G 

C PI TREE NOW IN ISAVE. 






MAIN6170 

671 

C 






MA1N618C 

572 

IF( IPRINT.LT. 3) SO TO 50 QC 






MAIN619Q 

673 

CALL TRPRNT( IrlSAVEf JSlt 6 ) 






MAIN62CC 

674 

GO TO 5000 






MAIN6210 

675 

C 






MAIN622C 

576 

C GENERATE SENSE TREE. 






MAIN6230 

677 

C 






MA1N624C 

678- 

77CC ISUC-INJEP 






MAIN6250 

679 

CALL TREE($896) 






MAIN626G 

680 

C 






MAIN6270 

681 

C INITIALIZE 






MA1N628C 
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682 

C 



MAIN6290 

663 


DC 7710 


MAIN63CC 

S84 


NB0UT(-I)rC 


MAIN631C 

665 


NBNUMa )r0 


MA1N632C 

686 

7710 

CONTIMUE 


MAIN6330 

687 

C 



MAIN634C 

683 

C 

ELIMINATE PATHS IN SENSE TREE IF THE> CONTAIN THE STANDBY 

BLOCK. 

MAIN6350 

669 

C 



MAIN636C 

69C 


NBOUT( 1 )rlRB( IN0£P»2> 


MAIN6370 

SSI 


NEOTMX-1 


HAIN636C! 

692 


CALL PATH0T(S998f ISAVEf JSl) 


MAIN6390 

693 


IF (IPRINT.LT.3) 6C TO 7720 


MAIN640C 

594 


CALL TRPSNTd »ISAVE* JSlt 7 ) 


MAIN6410 

695 

C 



MAIN642C 

696 

C 

ZERO OUT SENSE aLOCK NUMBER IN SENSE TREE- 


MAIN6430 

697 

C 



MAIN644C 

698 

7720 

NaNUMdJrINOEP 


MAIN645Q 

699 


NBNKAX-1 


MAIN646C 

7C0 


CALL NUMCUT(S998f ISAVEf JSl) 


MAINS47a 

7C1 


CALL SFLEFT(i296 fISAVE. JSl) 


MAIN648C 

702 


CALL 0UP0UT(S393f ISAVE# JSD 


MAIN6490 

7C3 


IF(IFRINT.LT-3) CO TC 5CCC 


MAIN65CC 

704 


call TRPRNni flSAVEf JSlt 8 ) 


MAIN6510 

705 

C 



MAIN652C 

706 

C 

CALCULATE PI, THE PRCBASILITY YOU NEED SENSE TREE* STORE 

IN PSNS. 

MAIN653Q 

7C7 

C 



MAIN654C 

708 

5CCC 

T-G. 


MAIN6550 

7C9 


NC W=D 


KA1N656C 

71C 


K = 0 


MAIN6570 

711 

9CG0 

CONTINUE 


MAIN658C 

712 


GO TO < 51C0*52C0*530a# 5400# 5500) fIDIST 


MAIN659Q 

713 

c 



MAIN66GC 

714 

5100 

CO 5101 I-1*N6MAX 


MAIN6610 

715 


DPREC--T*TL(I) 


MAIN6B2G 

716 


PROBK I )=R0(1 ) *OEXP (DPREC) 


MAIN6630 

717 

5101 

CONTINUE 


MAIN664C 

718 

C 



MAIN6650 

719 

C 

KRE ARRAY CONTAINS STANDBY BLOCKS THAT HAVE HAD THEIR PROBABILITY 

MAIN666G 

720 

C 

PREVIOUSLY calculated. 


MAIN6670 

721 

C 



MAIN668C 

722 


DO 5110 J=1*KRS1 


MAIN&890 

723 


IF(KRfc( J).EO.C) GO TO 511C 


MAIN670G 

724 


NB=KR3( J) 


MAIN67XG 

725 


PRCBl(NB)=RT(NBfK-d) 


MAIN672C 

726 

5110 

CONTINUE 


MAIN673C 

727 

c 



MAIN674C 

728 

c 

SUBROUTINE SYSP CALCULATES THE SYSTEM PROBABILITY - RETURNS TOT- 

MAIN6750 

729 

c 



KA1N676C 

730 


CALL SYSP 


MAIN8770 

751 

c 



MAIN678C 

732 


IF(NOW.EQ.I) GO TO 9001 


MAIN6790 

733 


K = K+1 


MAIN68CC 

734 


TTOT*K)/NSTEP 


MAIN681C 

735 


FSNS ( K ) -1 .CD'^C - TOT 


MAIN6820 

736 

c * 



MAIN&830 

737 

c 

PREVENT DIVISION BY ZERO WHEN 16 SIGNIFICANT FIGURES LOST 

. 

MAIN684C 

738 

c 



MAIN6850 

739 


IF CPSNS(K ).6T. D.CD+C ) 60 TC 5160 


HAIN686C 

740 


IF( K-EO-1 ) GO TO 5150 


MAIN&870 

741 

c 



MAIN688C 

742 

c 

NEXT STATEMENT ELIMINATES DIVISION BY ZERO LATER WHEN CALCULATING 

MAIN689C 

743 

c 

RT. THE NUMBER USED RESULTS FROM SUBTRACTING .9999999999 

99899999 

MAIN69CD 

744 

c 

FROM i-CD+G JUST BEFORE LOSING ALL SIGNIFICANT FIGURES- 

THIS WILI 

LMAIN6910 

745 

c 

CAUSE AN INTL6RATICN ERROR IN CALCULATING RT* THE RELlAfcl 

LITY OF 

MAIN692C 

746 

c 

A SENSE BLOCK CR SENSE TREE. THE SYSTEM RELIABILITY MAY 

NOT BE 

MAiri6930 

747 

c 

EFFECTED IF IT APPROACHES l.C. 


MAIN694C 

743 

c 



MAIN6950 

748 


FSNS(K) - +.i734723475S768C7C9D-C17 


KAIN696C 

750 

516C 

IF( K.LE.NSTEP ) GO TO 5 1C C 


MAIN6970 

751 


ISTEP=NSTEP+1 


MAIN698C 

752 


IF( IPRINT. N£. 4 ) GO TO 6CC0 


MAIN6990 

753 


PRINT 518C* CPSNE (K ) *Kr 1, ISTEP ) 


MAIN7C0C 

754 

5190 

FCRMATC 5026.18 ) 


KAIN7C10 

755 


G C T C 6 C 0 0 


MAIN7020 

756 

5200 

CONTINUE 


MAIN7030 
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757 

758 

759 

760 

761 

762 

763 

764 

765 

766 

767 

768 

769 

770 

771 

772 

773 

774 

775 

776 

777 

778 

779 

780 

781 

782 

783 

784 

785 
788 

787 

788 

789 
79C 

791 

792 

793 

794 

795 

796 

797 

798 

799 
SCO 
801 
802 
303 
8C4 
805 
8C6 
307 
8C6 
309 
810 
811 
812 

813 

814 
315 
816 
317 
818 
319 
820 
821 
822 

823 

824 
325 
826 

827 

828 
'32 9 

630 

331 

832 


MAIN704C 
MAIN7050 
MAIN706C 
MAIN7070 
HA1N706C 
MAIN7C90 

numerically integrate TLE reliability of the standby blocks. MAIN710C 
TIME VARIABLE - JT = O-MISSION TIMEITTOT). KT = 1-101 INCREMENTS. MAIN7110 

MAIN712C 
MAIN7130 
MAIN714G 
MAIN715Q 
MAIN716C 
MAIN7170 
MAIN718C 
MAIN7190 
MAIN72CC 
MAIN7210 
MAIN7220 
MAIN723Q 
MAIN724C 
MAIN7250 
MAIN7260 
MAIN727C 
MAIN728C 
MAIN7290 
MAIN730C 
MAIN731Q 
MAIN732C 
MAIN733Q 
MAIN7340 
MAIN7350 
MAIN736C 


5300 CONTINUE 
5400 CONTINUE 
5500 CONTINUE 

C SEE COMMENT STATEMENT AFTER 4500. 

GO TO 6000 
C 
C 
C 
C 
C 

6QCC 60 TO (61GCf64D0> tISTDE V 
C 

C METHOD FOR ONE STANDBY ELOCK (ISTD5 Y= II . 

C 

610D No-IRBdNDEPfE > 

RT( N8»l )::RC(N&) 

GC TO (6iC2»62C2*63C2#64C2t65G2lfIDIST 


6102 ISTEP-NSTEP41 

00 6170 KT-2»ISTEP 

DFREC - (-TT0T*(KT-1)/NSTEP >*TL(NB1 

DPRECi:: (-TT OT* ( KT-1 1 /NSTEPI ♦TLSCINOEP) 

PINT- RC(NB) *PSNS(1I ♦DEXP ( DP RE C ) ♦ S WP RO B ( I NDE P ) • DEXP ( DPRE Cl I 
JTEND-KT-1 
DO 6150 JT-ltJTEND 

0PREC2 = C (-TT0T*(KT-1)/NST£P)» TL(NB) I ♦( (-TT OT* ( JT- 1 1 /NS TEP ) ♦ 

1 (TLD(N6)-TL(NBI I I 

0PREC3z( (-TT0T»(KT-1)/N5TEP)* TL(NB) I ♦( (-TTOT*( JT I /NSTEP) ♦ 

1 ( TLD<NE)-TL(NB) ) ) 

0PREC4-( (-TTCT*{KT-1)/NSTCP)*TlS CINOEP) ) ♦( (-T T OT • ( JT- 1 1 /NS TEP I ♦ 

1 (TLDSdNDEPl-TLStlNDEP n I 

0PR£C5 = ( (-TT0T*(KT-1)/NSTEP)*TlS dNOEPi ) ♦( (-TTOT^IJT I /NSTEP)* 

1 (TLDSdNDEPI-TLSdNDEP m 
YAVE=C ( 0£XP(DPR£C2)*0EXP <0PR£C4) )♦( DEXP I DP REC 3 1 • DEX P( DPREC 5 I • I /2DQMAIN7370 
PINT = PINT •* RC(NB)*(PSNS(JT41)-PSNS( JTM ♦ YAVE * SWPRCBdNDEP) MAIN738C 


6150 

5170 

5888 

5202 
6302 
6402 
6502 
C SEE 
C 
C 
C 
C 
C 

6400 


CONTINUE 

RT<NBfKT)rFINT/PSNSC KT) 

CONTINUE 

IF CIPRINT.NE.4 ) GO TO 6C1C 

PRINT 5838f <RT(N3#KT)t KTritlSTEP) 

FCRMAT(8E15.9) 

60 TO 601C 

CONTINUE 

CONTINUE 

CONTINUE 

CONTINUE 

COMMENT STATEMENT AFTER 4500. 


METHOD FOR MORE THAN CN E STANDBY BLOCK ( IS TDB Vr2 1 
RC FOR ALL STANDBY 8L 0 CK S ( 5 1- 65 ) = 1 

N3=NST03Y( NSBYMX) 

RT(NB»l)ri.O 

90 TO C61C3f62C3f6303t64C3t65C3) flDIST 


6103 ISTEP-NSTEP41 

DC 65C0 KT-2.ISTEP 

DPREC = STBYPRtLSAVEf LJSlfROrTLfRTfKT f OfTTOT# NSTEPtTLO) 

DPRECI- (-T TOT* CKT-1) /NSTEP )*TLS (INDEP) 

PI NT-PSNSC D* DPREC* SWPROBt INDEP) *OEXP (DPRECI) 

JTEND-KT-1 
DO 6450 JTrifJTEND 

OF RE C4=( (-TTCT* (KT-1 )/N STEP ) *TlS (INDEP ))♦( (-T TOT* (JT- 1 l/NSTEP)* 

1 (TLOS(INOEP)-TLS(INOEP) )) 

DPRE C5-( ( -TTCT* (KT-1 l/NSTEP ) *TLS (INDEP ) ) ♦( (-TTOT* ( JT l/NSTEP)* 

1 (TLDS(lNOtP)-TLS(INOEP) )) 

JjT= jT-1 

YAVE=( ST6YPR( LSAVEtLJSl# ROt TL r RT # KT » jT • TT OT *NSTE Pt TLD ) *OEX P (DPREC5M AIN7720 
1 ) 4 STBYPR (LSAVEfL JSl tR C t T L tR T t K T # J JT t T TOT tNSTEP tTLD) *DEXP (DPRECMA IN773 C 

24)1/2. CD40 MAIN7740 

PINT :: PINT 4 ( PSnS( JT4 1 ) -PS NS t JT ) ) ♦ YAVE * S « PR OB ( INDEP I MAIN7750 

6450 CONTINUE MAIN7760 

RT(NB.KT)=PINT/PSNS( KT) MAIN777C 

6500 CONTINUE MAIN7780 

IF(IPRINT.NE.4) GO TO 6C1C MAIN77SC 


MAIN7390 

MAIN74CC 

MAIN7410 

MAIN742C 

MAIN7430 

MAIN744C 

MA1N74SC 

MAIN746Q 

MA1N747C 

MAIN748Q 

MAIN749C 

MAIN75C0 

MAIN751C 

MAIN752C 

MAIN753C 

MAIN7540 

MAIN7550 

MAIN7560 

MA1N757C 

MAIN7580 

MAIN759C 

MAIN7600 

MAIN761C 

MAIN762Q 

MAIN7630 

MAIN7640 

MAIN7&5C 

MAIN7660 

MAIN767C 

MAIN7680 

MAIN769C 

MAIN77C0 

MAIN771G 
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833 


PRINT 5863» <RT(N3tKT)» KTri,IST£P) 


MAIN78C0 

634 


GO TO 6C1C 


MAIN781C 

835 

62C3 

CONTINUE 


MAIN7820 

636 

6303 

CONTINUE 


MAIN783C 

337 

6403 

CONTINUE 


MAIN7840 

636 

6503 

CONTINUE 


MAIN7850 

839 


GO TO 6010 


MAIN7860 

84C 

c 



MAIN787C 

341 

C 



MAIN7880 

642 

C 

KRe ARRAY CONTAINS STANCBY BLOCKS WHOSE RELIABILITY HAS BEEN 

MAIN789C 

843 

C 

CALCULATED. 


MAIN7900 

644 

C 



MAIN7910 

845 

8C1C 

00 6GG4 N::ltL5 


MAIN7920 

646 


IF (KRE(N).EQ.C) 60 TO 6005 


MAIN793D 

847 

6CG4 

CONTINUE 


MAIN794Q 

848 

6DC5 

KRB(N>-NE 


MAIN795C 

849 


KRSl^N 


MAIN796C 

850 

C 



MAIN797C 

851 

c 

CHECK IF ALL SENSE BLOCKS ACCOUNTED FOR* IF 

SO GO TO 8000. 

MA1N7380 

852 

c 



MAIN73SC 

853 


IQ::0 


MAIN800Q 

654 


DO 6006 I=i*Ml 


MA1N801C 

855 


IF( IQ.EU.I ) GO TO 6007 


MAIN8020 

656 


IF(IS(I).EQ.1NDEP) 10=1 


KAIN803C 

857 


GO TC 6Q16 


MAXN8040 

858 

6007 

ISCI‘1>=IS(I ) 


MAIN8050 

853 

6Q16 

CONTINUE 


MAIN8060 

660 


IF(ISa) .EQ.O) GO TO 6008 


MAIN807C 

881 

6CG6 

CONTINUE 


MAIN8080 

662 

60:0 6 

NSR=NSR-1 


NA1N8G9C 

863 


IFd.EQ.2) GO TO 8000 


MAIN81C0 

864 


JJJ= JJJ+1 


MAIN811C 

865 


IFC JJJ.GE.M4) GO TO 8000 


MAIN8120 

866 


60 TO 7CCC 


MA1N813C 

867 

C 



MAIN8140 

866 

c 

FINALIZE VARIABLES AND CALCULATE RELIABILITY 

OF TOTAL SYSTEM. 

MAIN815C 

869 

c 

SET NOW = 1. 


MAIN8160 

870 

c 



MAIN817C 

871 

8CCC 

T=TTOT 


HAIN8180 

872 


K=NSTEP 


MAIN819C 

873 


N0W=1 


MAIN8200 

874 

c 



MAIN821G 

875 

c 

DUPLICATE JSAVE INTO ISAVE. 


MAIN8220 

876 

c 



HAIN823G 

377 


CALL OUPTREdSAVEf JSl *J SA VE* I JS 1) 


MAIN8240 

876 


GO TC 3000 


MAIN825C 

879 

c 



MAIN826Q 

88C 

c 

PRINT RESULTS IF RO = 1* SUPPRESS PRINTING 

OF RC. 

MA3N827C 

861 

c 



MAIN8280 

682 

3001 

PRINT 806 


MAIN829C 

833 

306 

FCRMATdHl ) 


MAIN8300 

684 

c 



MAIN831C 

885 


IF( NS8YMX.LE.C) GO TO 700 


MAIN832C 

886 

c 



MAIN8330 

887 

c 

PRINT STANDBY BLOCK S ( 5 1- 65 J AND THEIR TREES. 


MAIN834Q 

688 

c 

STORE ALL DIFFERENT STANDBY ELOCKS IN THESE 

TREES IN NRB ARRAY. 

MAIN835C 

389 

c 



HAIN8360 

8S0 


END FILE 1C 


MAIN837G 

391 


REWIND 10 


HAIN838a 

692 


DC 530 N=1*NS&YMX 


MAIN833C 

893 


READ(IO) NS* LJSl 


MAIN8400 

894 


DC 5C0 I=1#M2 


KAIN841C 

895 


READ« 10) i LSAVEt J.D* J = 1 *M1) 


MAIN8420 

896 

SCO 

CONTINUE 


MAIN843C 

397 


IFdPRINT.LE. 1) 60 TO 52 5 


MAIN844C 

838 


PRINT 510* NS 


MAIN8450 

893 

5 1C 

FORMAK •OSTANOBY TREE REPLACED BY BL0CK»I3) 


MAIN8460 

900 


CALL TRPRNT(5fLSAVE*LJSlt 0 ) 


MAIN847C 

9C1 

525 

CALL TRENUMaSAVEfLJSl) 


HAIN8480 

9C2 


DC 520 K=l*NbNMAX 


MAIN849C 

303 


NRS1=NR31*-1 


MAIN850Q 

904 


NRB(NRB1)=NBNUM(K) 


MAIN851G 

905 

52C 

CONTINUE 


MAIN8520 

306 

530 

CONTINUE 


MAIN8530 

907 

C 



MAIN8540 

908 

C 

PRINT THE ORIGINAL TREE WITH STANDBY BLO CKS C 51-65 )/( L1«L4 ) . 

HAIN855G 
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9CS 

C 




SIC 



IFaPRINT.LE.l) GO TO 7 CC 

MAIN8565 

911 



PRINT 5o0f Llf L4 

MAIN857C 

912 


SoC 

FCRMAT( *10RI6INAL PROEAEILITY TREE WITH STANDBY BLOCKS (*I3t 

MAIN858C 

913 


1 

. * THRU»l3f* ).•) 

MAIN859Q 

S14 



CALL TRPRNT(5*I5AV£» JSlt C ^ 

MAIN86G0 

915 



PRINT 8C6 

MAIN8610 

S16 

C 



MAIN862C 

917 


700 

PRINT 701 

MAIN8630 

916 


701 

FCRKATI* ACTIVE F/R DORMANT F/R R-INITIAL 

MAIN8B4C 

919 


X 

. RELIABILITY* ) 

MAlNd&SO 

520 



00 7C4 I-ltL4 

MAIN868C 

921 



IF(I.L£*M1) 60 TO 706 

MAIN867C 

922 



IF tPRCBKI) .LE.C) 60 TO 704 

MAIN368C 

923 



PRINT 703» If PRGBKI) 

MAIN869C 

924 


7G3 

FORM AT ( • DL0CK*I3f44XfD17 .71 

MAIN87CC 

925 



GO TO 704 

MAIN871C 

926 

C 



MAIN872C 

927 

C 


IF IN NR3 ARRAYf HAS BEEN REPLACEOf THEREFORE NO RELIABILITY. 

MAIN873Q 

928 

c 



MAIN874C 

929 


706 

DO 702 J=lfNRBl 

MAIN875C 

S3C . 



IF (I.NE.NRB( J) ) GO TO 7 C2 

MAIN876C 

931 



IF( RC(I) .6E.1.0) GO TO 711 

MAIN8770 

932 



PRINT 7C7t If TL(I)f TLD(U» RG(I) 

MAIN8780 

933 



30 TO 7C4 

MAIN879G 

934 


711 

PRINT 709t If TL(I)fTLD(Il 

MAIN88CC 

335 



60 TO 704 

MAIN8810 

936 


702 

continue 

MAIN882C 

937 

c 



MAIN3830 

938 

c 


DC NOT PRINT SENSE BLOCKS SINCE TLtD-Cf TLDCI)=Gf AND R0(I)=1 

AREMAIN884C 

939 

c 


NOT RELATED TO THE SENSE SWITCH OPTION PARAMETERS AND ARE ONLY 

. MAIN8850 

940 

c 


USED TO SET SENSE BLOCK PROBABILITY - l.G IN TREE PATHS. 

MAIN886C 

941 

c 



MAIN8870 

942 



DO 715 J-lfNNSR 

KAIN886C 

943 



IF( II5( J) .£Q. I) GO TO 7C 4 

MAIN889Q 

944 


715 

CONTINUE 

MAIN890C 

945 

c 



MAIN891C 

946 



IF (PROEKI J.LE.O. ) GO TC 704 

MAIN892G 

947 



rr(RC(I)-l.)7C5f7u8»705 

MAIN3S3C 

948 


705 

PRINT 7C7f If TUDf TLDCDf RC II) t PROElCl) 

MAIN894C 

949 


7C7 

FORMAT!* BL0CK*I3f 2E16.7 f Fl 2 .7 t 017. 7 ) 

MAIN895G 

950 



GO TC 704 

MAIN896C 

951 


7ca 

PRINT 733f If TLlDf TUJIDf PROBKI) 

MAIN897C 

952 


709 

FORMAT!* BL0CK*I3f 2E16.7 fl2XfDl7.7) 

MAIN898C 

353 


704 

CONTINUE 

MAIN3390 

954 

c 



MAIN9CCC 

355 

c 


PRINT SENSE SWITCH OPTIONS. 

MAIN901C 

956 

c 



MAIN9C2C 

857 



IF! NNSR.EG.C) GO TO 740 

MAlN9C3a 

958 



PRINT 72C 

MAIN9C4C 

959 


72C 

FORMAT! ‘CSENSE SWITCH ACTIVE F/R DORMANT F/R PROBABILITY *) MAIN9C50 

960 



DC 73C JritNNSR 

MAIN9C6C 

361 



NS-IIS! J) 

MAIN9070 

962 



PRINT 725f NSfTLS(NS) fTLDSINS) f SWFRCBINS) 

MAIN9C8C 

963 


725 

FORMAT!* 9L0CK*I3f 2E16.7 f Fl 2 . 7 f E 1 6 . 7 ) 

MAIN9C9Q 

964 


730 

CONTINUE 

MAIN91CC 

965 


740 

RzTCT 

MAIN9110 

966 



K-“LC616!i.-R) 

MAIN912C 

967 



K=K+NSIG 

MAIN9130 

968 



IFIK.GT.6) K-6 

MAIN9131 

969 



V!12)::VK(Xf‘2) 

MAIN9132 

970 



V! 14 )=VK ! K) 

MAIN9135 

371 



WRlTE(bfV) TTOTfR 

MAIN9137 

972 



PRINT 6Co 

MAIN914C 

973 

c 



MAIN9150 

974 

c 


MODIFY PROBABILITY DATA IF LAST = 7f GO TC 9£0. 

MAIN916C 

375 

c 



MAIN3170 

976 



IF !L AST.LT.7 ) 60 TO 998 

MAIN918C 

977 



IF! LAST. EO. 7) GO TO 38C 

MAIN9190 

978 



IF !L AST.EQ.8 ) GO TO 939 

MAIN92CG 

973 



IF! LAST. EO. 9) GO TO 1C6S 

MAIN9210 

980 



60 TC 996 

MAIN922C 

931 

c 



MAIN9230 

S62 

c 


INITIALIZE variables FOR RECALCULATION WITH NEW PARAMETERS. 

MAIN924C 

983 

c 



MAIN9250 

984 


S6C 

DO 99C I^lfMi 

MAIN526C 
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985 

990 

I3( DzllSi I) 

S66 


N5R-NNSR 

987 


I3UC= JSUC 

988 


GO TO 997 

989 

C 


S5C 

C 

ERROR RETURN FROM SU cRO L'T INE S . 

931 

c 


992 

596 

READ 2707 . LAST 

993 

2707 

FORMAT (73XrIl) 

594 


IF(LAST-8) 9S8t953.1C&6 

995 

1058 

END 


MAIN9270 

MAIN928C 

MAIN3290 

MAIN931G 

MAIN932Q 

READ TO NEW DIAGRAM DATA DECK. MAIN933C 

MAIN934G 

MAIN935C 

MAIN9360 

MAIN937C 

MAIN938Q 
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IREE»F TREE.SFLEFT 


1 

2 

C 

c 


SUBROUTINE SFLEFI SHIFTS PATHS IN ISAVE TO THE'lEFT' 
NB AT LEFT-HOST. FILL IN ZEROS AFTER SHIFTIN3. IF 

Until nonz ero 
50 ZEROS PUTS 

slft 

SLFT 

10 

20 

3 

c 


P«TH NUM3ER IN NPUP ftNO ELIMINATES IF BF CALLING PATHUP. 

slft 

30 

4 

c 


yRiTTEN BY RON ECKSTEIN Wl/70 


_iLFr 

40 

5 

c 




slft 

50 

6 



SUBROUTINE SFLEFT ( I a IS A VE s JS 11 


_slft__ 

SO 

7 

c 




slft 

70 

3 



COMMON/ ALLSUE/M lfM2»M3iM4tM5fM&M7.LltL3.L4*L5#L6 


_SL.FJ__ 

80 

9 



COMMON/SMSU3/NPUP ( 200) tNP UP MX •N30UT 1 501 • N3 OT MX . N 3IN t 50 ) # N8 INMX • 

slft 

90 

10 


1 

N3NUMI 50) fNoNM AX 


slft 

100 

1 1 



DIMENSION ISAVECHl. M2» 


slft 

110 

12 

c 




SUJ 

!^D_ 

1 3 

c 


initialize 


slft 

1 30 

14 



DO 10 I=lfM2 


.SL.FJ_. 

J4_0_ 

15 



NPUPI I )=0 


slft 

150 

16 


10 

CONT INUE 


SLFI. 

J6.Q_ 

1 7 



L=0 


slft 

1 70 

IS 

c 




SLJT 

130 

19 



00 ICO IrifJSl 


SLFT 

1 90 

20 



ILEFTrO 


..SLF_T_. 

_2Q.O. 

21 



00 90 J=liMl 


slft 

210 

22 



IFITSAVEC J»I l-NE.D) GO TO 70 


SL-FJt 

_22Q 

23 



ileft=ileft»i 


slft 

2 30 

24 



GO to 90 


^SL.E-T„ 

.2.413. 

25 


70 

IF( IL£FT,£Q,C ) GO TO 90 


slft 

2 5n 

26 

c 




_SLFT„ 

.Z6a 

27 

c 


MOVE N3 TO LEFT 


slft 

2 70 

28 

f 




SLFI,. 

28.0 

29 



M=J-IL£FT 


slft 

2 90 

30 



ISAVE (Mfl )=l SAV£( J>I 1 


SLFT. 

300. 

31 


90 

CONTINUE 


slft 

310 

32 



IF CILEF T.NE.MI 1 GO TO 1 20 


SLFT. 

320 

33 



L::L+1 


slft 

3 30 

34 



NPUPCL)=I 


slft 

340 

35 



30 TO 100 


slft 

3 50 

36 

F* 

w 




-SLFT 

360 

37 

c 


ZERO OUT REDUNDANT 3L0CK NUMBERS. 


slft 

3 70 

39 

.. C. 




slft 

380 

3 9 ‘ 


120 

I9EG = MI-IL£FT M 


SLFT 

3 90 

40 



DO 150 K=I3EG»MI 


SLFT 

400 

4 1 



ISAVE(K#I)=0 


slft 

4 10 

42 


ISO 

CONT INUE 


SLFT 

420 

V3 " 


100 

CONTINUE 


slft 

4 30 

44 



NPUPMXZL 


SLFT 

440 

1,5 



IF( NPUPMX.EO. 0) GO TO 1000 


slft 

4 50 

46 



CALL PATHUP ($1100. ISAVE fjSl) 


SLFT 

460 

ei7 


1000 

RETURN 


slft 

"4 70 

49 


lino 

PRINT llO! 


slft 

480 

49 


i ini 

FORMAT!" ERROR IN PAThIP CALLED FROM SFlEFT,") 


slft 

4 30 

50 



RETURN 1 


._s_l f_i 

500 

51 



END 


SLFT 

510 


16 
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^♦FTREE 

.OUPOUT 

1 

C 

2 

C 

3 

C 

4 

C 

5 


6 

C 

7 


8 


9 


10 


11 

C 

12 

C 

13 

C 

14 


15 


16 

10 

17 


18 

C 

19 


20 

c 

21 

c 

22 

c 

23 


24 


25 


26 

85 

27 

c 

28 

C 

29 

c 

30 

90 

31 

95 

32 

C 

33 

C 

34 

c 

35 


36 


37 


36 

80 

39 

75 

40 


41 

70 

42 

c 

43 

c 

44 

C 

45 


46 


47 

60 

48 


49 

100 

50 


51 


52 


53 

1000 

54 

HOC 

55 

1101 

56 


57 



SUBROUTINE OUFOUT STORES PATH NUMBERS THAT ARE DUPLICATES IN NPUP DPOT 


AND ELIMINATES THEM BY CALLING PATHUP* DPOT 

WRITTEN BY RON ECKSTEIN 1/1/7C DPOT 

DPOT 

SUBROUTINE DUP 0 UT ( $ # ISA VE t JS 1 ) DPOT 

DPOT 

CCMM0N/ALLSUB/MltM2fK3f M4 tM5 tM 6 f M 7 f LI t L3 t L4 t L5 *L6 DPOT 

COM MON/ 5 MS UB/ NPUP ( 200) tNPUP MX t NBOUT (50) »NB OTMX t NBIN ( 50 ) » NBINHX • DPOT 

1 NENUM(5C ) tNBNMAX DPOT 

DIMENSION ISAVE(M1# M2) OPOT 

DPOT 

INITIALIZE DPOT 

DPOT 

DO 10 I=1#M2 DPOT 

NPUP (I )=0 DPOT 

CONTINUE DPOT 

LzC DPOT 

DPOT 

DO IDD I::ifJSl DPOT 

OPOT 

CHECK IF PATH ALREADY A DUPLICATE OF PREVIOUS PATH. DPOT 

OPOT 

IFCL.EQ.C) 6C TO 9C DPOT 

DO 85 J-1#L DPOT 

IF(NPUP(J).EG.I) GO TO lOD DPOT 

CONTINUE DPOT 

DPOT 

CHECK PATH AGAINST ALL HIGHER PATH NUMBERS. OPOT 

DPOT 

DPOT 

IF (II.GT. JSl ) GO TO IDD DPOT 

OPOT 

CHECK IF PATH ALREADY A DUPLICATE OF PREVIOUS PATH. OPOT 

OPOT 

IF (L.EG.O) GO TO 75 DPOT 

DO 80 J=i#L OPOT 

IF (NPUP( JI.EO.II) 60 TO 6C OPOT 

CONTINUE DPOT 

DO 70 K^ltMl DPOT 

IF( ISAVE(K#I) .NE.ISAVE (K f ID) GO TO 60 OPOT 

CONTINUE DPOT 

OPOT 

STORE DUPLICATE PATH NUMBER IN NPUP. DPOT 

OPOT 

L-L+1 DPOT 

NPUP(L):=II OPOT 

Iini + l DPOT 

60 TO 95 OPOT 

CONTINUE OPOT 

NPUPMX=L DPOT 

IF (NPUPMX.EG.O) 60 TO ILOC DPOT 

CALL FATHUPlS110D»lSAV£f JSl) OPOT 

RETURN DPOT 

PRINT liCi OPOT 

FORMAT(* ERROR IN PATHUP CALLED FROM DUPOUT.M DPOT 

RETURN 1 OPOT 

END DPOT 


1C 
20 
3C 
40 
5D 
60 
7C 
80 
90 
ICO 
110 
120 
130 
140 
15C 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
2 80 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
39 0 
4 GO 
410 
420 
430 
440 
450 
460 
470 
480 
49 C 
500 
51C 
520 
530 
540 
550 
560 
570 
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TREA^FTREE.PATHUF 


1 

C 


subroutine pathup takes path numbers to be 

ELIMINATED FROM NPUP 

PAUP 

10 

2 

C 


ARRAYf ELIMINATES THESE PATHS FROM ISAVEt 

MOVING UP REMAINING PATHPAUP 

20 

3 

c 


AS IT DOSSf AND ZEROS OUT REDUNDANT PATHS 

AT THE END OF THE TREE. 

PAUP 

30 

4 

c 


WRITTEN EY RON ECKSTEIN 1/1/70 


PAUP 

4C 

5 

c 




PAUP 

50 

6 



SUBROUTINE P ATHUP { $ 9 ISA \iE » JS 1 ) 


PAUP 

60 

7 

c 




PAUP 

70 

8 



CC MM ON/ ALLSUL/Mli M2 9 H3 9 K4#M5 9 M 6 9 M 7 tL 1 9 L3 t L4 # L5 9 L6 

PAUP 

80 

9 



COM MON/SMS U3/NPUP (200) tNPUP MX »N60UT ( 50) 9 N 3 OTMX t N BIN ( 50 ) . NB INMX # 

PAUP 

90 

10 



1 NENUM(5G) 9NBNMAX 


PAUP 

ICO 

11 



DIMENSION ISAV£(Ml9 M2) 


PAUP 

lie 

12 

c 




PAUP 

12 c 

13 

c 


CHECK IF THERE ARE PATHS TO 3£ ELIMINATED. 


PAUP 

13C 

14 

c 




PAUP 

14C 

15 



IF(NPUPMX.EQ-Q) GO TO 200 


PAUP 

150 

16 

c 




PAUP 

160 

17 



IUP=u 


PAUP 

170 

18 



DO IOC I=l9JSl 


PAUP 

180 

19 



DO 90 J-ltNPUPMX 


PAUP 

130 

20 



IF(I.NE.NFUP(J) ) GO TO SC 


PAUP 

200 

21 



IUP=IUP+1 


PAUP 

210 

22 



IF (lUP.GT.NPUPMX) GO TO llD 


PAUP 

220 

23 



GO TO IQC 


PAUP 

2 30 

24 


90 

CONTINUE 


PAUP 

24C 

25 



IF( lUP.EQ.O) GO TO 100 


PAUP 

2 50 

26 

c 




PAUP 

260 

27 

c 


MOVE PATH NUMBER H UP A NUMBER r lUP. 


PAUP 

2 70 

28 

c 




PAUP 

280 

29 


60 

M=I-IUP 


PAUP 

2 90 

30 



DO 40 K-I 9 MI 


PAUP 

30 C 

31 



ISAV£CKfM)rISAVE(K#I) 


PAUP 

310 

32 


40 

CONTINUE 


PAUP 

32 C 

33 


100 

CONTINUE 


PAUP 

3 30 

34 



IFCIUP.EQ.NPUPMX) GO TO 150 


PAUP 

340 

35 


lie 

PRINT 120# lUPfNPUPMX 


PAUP 

350 

36 


120 

FORMAT( »1ERRCk IN PATHUF AT IIC. IUP=»I4# 

• NPUPMX=»I4) 

PAUP 

360 

37 



RETURN 1 


PAUP 

370 

38 


15C 

IEND= JSl 


PAUP 

380 

39 

c 




PAUP 

3 9G 

40 

c 


RESET JSl 


PAUP 

400 

41 

c 




PAUP 

410 

42 



JS1= JSl-IUP 


PAUP 

420 

43 

c 




PAUP 

430 

44 

c 


ZERO OUT REDUNDANT PATHS. 


PAUP 

440 

45 

c 




PAUP 

4 50 

46 



ISTART-JSl+1 


PAUP 

460 

47 



DO 16D IziSTART# lENO 


PAUP 

470 

48 



CO 160 J=l9Ml 


PAUP 

480 

49 



I5AVE( JfI)=C 


PAUP 

4 90 

50 


16C 

CONTINUE 


PAUP 

5CC 

51 


200 

RETURN 


PAUP 

510 

52 



END 


PAUP 

520 


18 


JPL Technical Memorandum 33-513 



TREA»FIRE£.PATHOT 


1 

C 


subroutine: pathot compares numbers in nbout array with each 

NUMBERPAOT 

10 

2 

C 


IN each success path in isave and generates npup array with 

PATH 

PACT 

20 

3 

C 


NUMBERS to BE ELIMINATED IE ANY 

ABSOLUTE VALUES ARE EQUAL. 


PAOT 

30 

4 

C 


WRITTEN BY RON ECKSTEIN 1/1/70 



PACT 

40 

5 

C 





PAOT 

50 

6 



SUBROUTINE PA THO T ( $ t IS AV E » JSl ) 



PACT 

60 

7 

C 





PACT 

70 

S 



C0MMCN/ALLSUB/Ml*M2»M5»M4f M5f M6 

•M7fLlfL3fL4#L5#L6 


PAOT 

80 

3 



C0MM0N/SMSUB/NPUF (20C>iliPUPMXtNBCUTC5D) f NBCTMXtNBlN(5C>»NEINMX» 

PAOT 

90 

10 



1 NBNUMC5C) »NBNMAX 



PAOT 

100 

11 



DIMENSION ISAVE(Mlf M2) 



PAOT 

lie 

12 

C 





PAOT 

120 

13 

C 


INITIALIZE 



PAOT 

13C 

14 

C 





PAOT 

140 

15 



DO 1C I-lfM2 



PAOT 

15C 

16 



NPUPC I)rC 



PAOT 

160 

17 


1C 

CONTINUE 



PAOT 

17C 

18 



L-0 



PAOT 

180 

13 

C 





PAOT 

190 

2C 



DC IOC I-itJSl 



PAOT 

200 

21 



DO 90 J-1»M1 



PACT 

210 

22 



DO 80 K-lfNBOTMX 



PAOT 

220 

23 



IF(NeCUT(K).NE.ABS(ISAVECJtin)60 TO 6C 


PAOT 

23C 

24 



L = L+1 



PAOT 

240 

25 



NPUP<L)=I 



PAOT 

25C 

26 



GO TO IOC 



PACT 

260 

27 


ec 

CONTINUE 



PAOT 

27C 

23 


90 

CONTINUE 



PAOT 

280 

29 


ICC 

CONTINUE 



PAOT 

290 

3C 



NPUPHX-L 



PAOT 

300 

31 



IF (NPUPMX.EQ.C) GO TO ICCO 



PAOT 

310 

32 



CALL PATHUP(Sliaa#ISAVE» USD 



PAOT 

320 

33 

ICCC 

RETURN 



PAOT 

330 

34 

1100 

PRINT ilCl 



PAOT 

340 

35 

1101 

rORMATC ERROR IN PATHUF CALLED 

FROM PATHOT. M 


PAOT 

350 

36 



RETURN 1 



PAOT 

360 

37 



END 



PACT 

370 
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TREA*FTR£E.NUMOUT 


1 

C 


SUBROUTINE NUMOUT ZEROES OUT ALL NUMBER BLOCKS IN 

ISAVE THAT ARE 

NMOT 

10 

2 

C 


EQUAL IN ABSOLUTE VALUE TO ANY OF THE NBNUM ARRAY 

NUMBER BLOCKS. 

NMOT 

20 

3 

C 


IF 50 ZEROESf STORE PATH NUMBER IN NPUP AND E LIMINATE f PATHUP) - 

NMOT 

30 

4 

C 


WRITTEN BY RON ECKSTEIN 1/1/70 


NMOT 

40 

5 

c 




NMOT 

50 

6 



SUBROUTINE NUMOUU $ t ISA VE t JS 1 I 


NMOT 

&C 

7 

c 




NMOT 

70 

8 



COMMON/ ALLSUB/Ml*M2t M3 rM4 tM5 #M 6 t M 7 tLlt L3 » L4 t L5 tL6 


NMOT 

80 

9 



CCMM0N/SMSU3/NPUP* 200) tN PUP MX » N80UT < 50) » NB OTMX * NBIN (50 ) t N8 INMX # 

NMOT 

90 

ID 



1 NBNUM(50) tNBNMAX 


NMCT 

100 

11 



DIMENSION ISAV£(M1» M2) 


NMOT 

110 

12 

c 




NMOT 

12 0 

13 

c 


INITIALIZE 


NMOT 

130 

14 

c 




NMOT 

14 0 

15 



DC 10 I^lfM2 


NMOT 

150 

16 



NPUP (I )=D 


NMOT 

160 

17 


1C 

CONTINUE 


NMOT 

170 

18 



NPUPMX=:0 


NMOT 

18C 

19 



L = 0 


NMOT 

190 

20 

c 




NMOT 

200 

21 



DO ICC I^lfJSl 


NMOT 

210 

22 



IC0UNT=0 


NMOT 

220 

23 



DO 90 J^lfMl 


NMOT 

230 

24 



DO 80 K=lfNBNMAX 


NMOT 

240 

25 



IF( ISAVE( J fl) .EQ.O) GO TO 75 


NMOT 

2 50 

28 



IF (NBNUM( K) • EQ» AbS( ISAVEC Jtl )n60 TO 70 


NMOT 

260 

27 


30 

CONTINUE 


NMOT 

270 

28 



GO TO 90 


NMOT 

280 

29 


70 

ISAVEi J»I)::a 


NMOT 

290 

30 


75 

IC0UNT::IC0UNT+1 


NMOT 

300 

31 


30 

CONTINUE 


NMOT 

310 

32 

c 




NMOT 

320 

33 

c 


IF 50 ZEROES STORE PATH NUMBER IN NPUP. 


NMOT 

330 

34 

c 




NMOT 

34 0 

35 



IF( ICOUNT.lt. M l) GO TO 100 


NMOT 

350 

38 



L = L+1 


NMOT 

360 

37 



NPUPIL)::! 


NMOT 

370 

38 


ICO 

CONTINUE 


NMOT 

38 0 

39 



NPUPMX^L 


NMOT 

3 30 

40 



IF (NPUPMX.EQ.D) GO TO ICDO 


NMOT 

40G 

41 



CALL PATHUP($110G#ISAVEf JSl) 


NMOT 

410 

42 

IGCO 

RETURN 


NMOT 

420 

43 

1100 

PRINT 1101 


NMOT 

430 

44 

1101 

FORMAT!* ERROR IN PATHUP CALLED FROM NUMOUT.*! 


NMCT 

440 

45 



RETURN 1 


NMOT 

4 50 

46 



END 


NMOT 

46C 


20 
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.NEGOUT 

C 

C 

C 

c 


SUBROUTINE NEGOUT STORES FAILURE PATH NUMBERS ENDING WITH A 
NEGATIVE SIGN IN NPUP ARRAY# THEN CALLS PATHUP TO ELIMINATE- 
WRITTEN BY RON ECKSTEIN 1/1/70 

SUBROUTINE N EGO UT C $ » ISA VE • JS 1 I 


CCMM0N/ALLSUB/MltM2 #M3#KH # M5 #M6 tM 7 #L1 # L3 # L4 » L5 1 L6 
CCMM0N/SMSU3/NPUP (200) »NPUPMX »NBOUT (50> # NB OTMX # NBIN(50 1 # NBINMX # 
1 NBNUM(5C ) tNBNMAX 
DIMENSION ISAVE(M1# M2) 

INITIALIZE 

00 IG I::i#M2 
NPUP(I)=0 
: CONTINUE 
Lro 

DO ICC I = 1#JS1 
DO 8C J^itMl 

IF (ISAVEC Jfl )-EQ-C) GO TO 90 
3 CONTINUE 
J:iLl 

CHECK LAST NONZERO BLOCK NUMBER IN PATH FOR SIGN- 


IF(J.LE.l) GO TO 92 

IFC ISAV£( J-1# I) ) 91#92#X00 

L = L-H 

NPUPCD-I 

60 TO 100 

PRINT 95 

FORMAK *1ERRCR IN SUBROUTINE NEGOUT AT 92-M 
IF{ U.GT.l) RETURN 1 
PRINT 94 

FCRMATC* ISAVEdf 1)=C* ) 

GO TO lOOC 

CONTINUE 

NPUPMX-L 

IF( NPUPMX-EQ- 0) GO TO 1000 
CALL FATHUPC$11DG#ISAVE fJSl) 

RETURN 
PRINT 1101 

FORMAT! • ERROR IN PATHUP CALLED FROM NEGOUT. 

RETURN 1 

END 


NGCT 10 
NGOT 20 
NGCT 3C 
NGOT 40 
NGOT 50 
NGOT 60 
NGOT 70 
NGOT 80 
NGOT 90 
NGOT 100 
NGOT 110 
NGOT 120 
NGCT 130 
NGOT 140 
NGCT 150 
NGOT 160 
NGCT 170 
NGOT 18C 
NGCT 190 
NGOT 2C0 
NGOT 210 
NGOT 220 
NGCT 230 
NGOT 240 
NGOT 25C 
NGOT 260 
NGOT 270 
NGOT 28C 
NGOT 290 
NGOT 300 
NGOT 310 
NGOT 320 
NGOT 330 
NGOT 340 
NGCT 350 
NGOT 360 
NGOT 370 
NGOT 380 
NGOT 390 
NGOT 400 
NGCT 410 
NGOT 420 
NGOT 430 
NGOT 440 
NGOT 45C 
NGOT 460 
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TftEA»FTREE;.TRINTR 


1 

C 


SUBROUTINE TRINTR PUTS THE STANDBY BL0CK(51-651 IN PLACE OF 

THE 

TRTR 

10 

2 

C 


FIRST PATH OF THE STANDBY TREE IMBEDDED IN THE PATH OF THE 


TRTR 

2C 

3 

C 


ORIGINAL TREE IN JSAVE» THEN ELIMINATES ORIGINAL TREE PATHS 

THAT 

TRTR 

30 

** 

C 


CONTAIN 2ND-JS1 STANDBY PATHS. THEN PUTS THE NEGATIVE STANDBY 

TRTR 

4C 

5 

C 


6L0CK( 51-65)/ (L1-L4) IN AND ELIMINATES CORRESPONDING PATHS. 


TRTR 

50 

6 

c 


WRITTEN BY RON ECKSTEIN 1/1/70 


TRTR 

60 

7 

c 




TRTR 

70 

8 



SUBROUTINE T RIN TR ( $ • JSA VE tiJ SI # IS A VE t JSl t NST DB Y tNSBY MX tNSS ) 


TRTR 

80 

9 

c 




TRTR 

90 

10 



CCMMCN/ALLSUB/M1.M2 thSt M4 #M5 tM6 *M 7 #L1 # L3 * L4 t L5 *L6 


TRTR 

IDO 

11 



COMMON/SMSUB/NPUP (200) tNPUP MX t NBOUT (50) f NB OTMX • N BIN ( 50 ) * NBI NMX • 

TRTR 

110 

12 



. NBNUM( 50 ) f NBNMAX 


TRTR 

120 

13 



DIMENSION JSAVE(Mlf M2 ) r IS A VE ( Ml t M2) * NSTOB Y ( M4 ) 


TRTR 

130 

14 

c 




TRTR 

14 0 

15 

c 


INITIALIZE 


TRTR 

150 

18 

c 




TRTR 

16 0 

17 



JJ=C 


TRTR 

170 

16 

c 




TRTR 

18 0 

19 


5 

00 10 J=1»M2 


TRTR 

190 

20 



NPUP( J)=G 


TRTR 

20 0 

21 


10 

CONTINUE 


TRTR 

210 

22 



L = 1 


TRTR 

220 

23 

c 




TRTR 

2 30 

24 



DO 1000 I=1#IJS1 


TRTR 

24 0 

25 



IF(L.GT.JSl) GO TO 1010 


TRTR 

2 50 

26 



J=1 


TRTR 

260 

27 



K=1 


TRTR 

270 

28 



ISTART::C 


TRTR 

280 

29 


30 

IF(J.GT.Mi) GC TO 1000 


TRTR 

2 90 

30 



IF (JSAVE( J»I J.EQ.O) GO TC ICCC 


TRTR 

30 0 

31 



IF( JSAV£( Jf I) .NE.ISAVE (K #L ) ) GO TO 80 


TRTR 

310 

32 



IF (ISTART.NE.C) GO TC 50 


TRTR 

320 

33 



ISTARTrJ 


TRTR 

330 

34 



NPATH=I 


TRTR 

340 

35 


50 

KrK+1 


TRTR 

350 

36 



IF (K.6T.M1) GO TO IGO 


TRTR 

360 

37 



IFdSAVE(KfL) .Ea.O) 60 TO lOQ 


TRTR 

370 

38 


80 

J= J-H 


TRTR 

380 

39 



GO TO 33 


TRTR 

3 90 

40 


100 

IFIL.GT.I) GC TO 130 


TRTR 

40 G 

41 



IST0P-I5TART+K-2 


TRTR 

410 

42 



GO TC 300 


TRTR 

42 0 

43 

c 




TRTR 

4 30 

44 

c 


STORE original TREE PATFS THAT CONTAIN 2ND-JS1 STANDBY PATHS 

• 

TRTR 

44C 

45 

c 




TRTR 

450 

46 


130 

DO 150 J=1*M2 


TRTR 

460 

47 



IFI NPUP( J) .NE .0) GO TO 150 


TRTR 

470 

48 



NPUP( J)=I 


TRTR 

48C 

49 



NPUPMX-J 


TRTR 

490 

50 



GO TO 900 


TRTR 

50 C 

51 


150 

CONTINUE 


TRTR 

510 

52 



PRINT 160 


TRTR 

520 

53 


1 60 

FORMAT( dERROR IN TRI-TR AT 160. -PUP IS FULL.*) 


TRTR 

530 

54 



RETURN 1 


TRTR 

54 0 

5 5 


300 

IF(L.Ea.l) GC TO 304 


TRTR 

5 50 

56 



PRINT 301» L 


TRTR 

56C 

57 


301 

FORMAT( *1ERR0R IN TRINTR AT 300. L=*I4) 


TRTR 

570 

58 



RETURN 1 


TRTR 

580 

59 

c 




TRTR 

5 90 

60 

c 


IF NSS EQUAL ZEROt TRINTR IS BEING CALLED FROM STDBY2. OTHERWISE 

tTRTR 

60 0 

61 

c 


TRINTR IS BEING CALLED FROM MAIN PROGRAM TO REPLACE STANDBY 

TREES 

TRTR 

610 

62 

c 


PREVIOUSLY ASSIGNED A STANDBY BLOCK (51-65 )/( L1-L4 ) . 


TRTR 

620 

63 

c 




TRTR 

6 30 

84 


304 

IF (NSS.EQ.O) GO TO 305 


TRTR 

640 

65 



JSAV£(ISTART#NPATH)=NSS 


TRTR 

6 50 

66 



GO TO 325 


TRTR 

660 

67 

c 




TRTR 

670 

68 

c 


ASSIGN standby BLOCK NUMBER ( 51-65 ) ONLY ONCE 


TRTR 

680 

69 

c 




TRTR 

6 90 

70 


305 

IF(JJ.NE.C) GO TO 313 


TRTR 

70 0 

71 



DO 310 J=1»M4 


TRTR 

710 

72 



IF (NSTOEY( J) .EQ.O 1 GO TC 312 


TRTR 

720 

73 


310 

CONTINUE 


TRTR 

730 

74 



PRINT 315 


TRTR 

74 0 
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75 

76 

315 

77 

76 

312 

79 

C 

80 

C 

81 

C 

82 

C 

33 

84 

313 

35 

325 

86 

C 

87 

C 

88 

C 

89 

9C 

C 

91 

S2 

320 

93 

S4 

330 

95 

SCO 

96 

IDCO 

97 

C 

98 

C 

99 

IOC 

c 

101 

1C2 

1010 

1C3 

c 

104 

c 

105 

c 

106 

c 

107 

c 

1C8 

c 

103 

c 

IID 

c 

111 

c 

112 

1111 

113 

c 

114 

c 

115 

116 

117 

118 

c 

119 

450 

12G 

C 

121 

122 

123 

124 

125 

460 

126 

127 

500 

128 

C 

129 

c 

130 

c 

131 

SIC 

132 

c 

133 

c 

134 

135 

136 

137 

138 

139 

c 

140 

141 

520 

142 

525 

143 

144 

530 

145 

oCC 

146 

147 

2000 

148 

c 

149 

c 


FCRMATC ’iERROR IN TRINTR AT 315. NSTOBY IS FULL.M 
RETURN 1 

JJ-J 

NSBYKX=JJ 

PUT NEXT standby BLOCKS 51*65) IN PLACE OF THE FIRST PATH OF THE 
STANDBY TREE IMBEDDED IN THE PATH OF THE ORIGINAL TREE. 

NSTD8Y( 

JSAVE CIST ART fNPATH) = NSl DB Y (J J) 

ISTART-ISTART+1 


TRTR 750 
TRTR 7SC 
TRTR 77G 
TRTR 78C 
TRTR 790 
TRTR 800 
TRTR 81Q 
TRTR 820 
TRTR 83C 
TRTR 840 
TRTR 850 
TRTR 86 C 


ZERO OUT 2NCHNTH NUMBER BLOCKS IN FIRST PATH IN ISAVE IMBEDDED IN TRTR 870 


JSAVE. 

IF (ISTART.GT.ISTOP) 60 TO SOD 
DO 330 M-iSTART»ISTOP 
JSAVE(M.NPATH)=Q 
CONTINUE 

CALL SFLEFT($110Dt JSAVE tiJSl ) 

L-L+1 

CONTINUE 

CHECK IF NO MORE STANDBY TREES IN ORIGINAL TREE. 

IF (ISTART.EQ.O) 60 TO 1111 
CALL PATHUPf SilOOf JSAVEf USD 
60 TO 5 


TRTR BSC 
TRTR 890 
TRTR 300 
TRTR 910 
TRTR 920 
TRTR 930 
TRTR 940 
TRTR 950 
TRTR 960 
TRTR 970 
TRTR 980 
TRTR 990 
TRTRIOOO 
TRTRlOlO 
TRTR1020 
TRTR1030 
TRTR1C4C 
TRTR1050 
TR TR106C 


PUT NEGATIVE STANDBY 3L0 CK ( 51-65 )/ C L1-L4) IN PLACE OF FIRST NB 
THAT IS equal IN ABSOLUTE VALUE TO ANY OF THE NBNUM ARRAY OF THE 
TREE THAT THE STANDBY BLOCK (51-651/ CL1-L4) REPLACES. ZERO OUT THETRTR1070 
FOLLOWING NB IF THEY ARE EQUAL IN ABSOLUTE VALUE TO ANY OF THE TRTR1C8G 


NBNUM ARRAY. CALL SFlEFT TO ELIMINATE AND RESTORE PATHS. 


CALL TRENUMCISAVE. JSl J 

ASSIGN STANDBY BLOCK ( 51 -6 5 )/ (L 1-L4 ) ONLY ONCE. 

IF(NSS.EG.C) GO TO 450 
NS=NSS 
60 TO 460 
NS=NSTDSY( JJ) 

00 20C0 I-1#1JS1 

DO 600 J=l#Mi 

00 SOD K-1»NBNMAX 

IF ( JSAVEIJfl ).EQ.C) 60 TO 2CCD 

IF( N3NUH( K ) .EG. AOS ( JSAVE ( Jf I) ) ) GO TO 510 

CONTINUE 

GO TO 6GC 

REPLACE WITH NEGATIVE STAND3T BLOCK ( 51-65) / (L1-L4) . 
JSAVE( JtI)--NS 

ZERO OUT ANY NB FOLLOWING IF EQUAL TO NBNUM ARRAY. 

LSTART^J^l 

DO 530 M-LSTARTfMl 

DC 520 N::1»NBNMAX 

IF ( JSAVKMtl ).EQ.C) GO TO 20GC 

IF( NBNUfM N) .EQ.ABSC JSAVE (Mf I) ) ) GO TO 525 

CONTINUE 

60 TO 20CD 

JSAVE ( Mtl )-0 

CONTINUE 

GO TO 2000 

CONTINUE 

CONTINUE 

CALL SFLEFTtSllOCf JSAVEt USD 
ELIMINATE DUPLICATE PATHS. 


TRTR1090 

TRTRllOC 

TRTRlllO 

TRTR112G 

TRTR1130 

TRTR114C 

TRTR1150 

TRTR1160 

TRTR1170 

TRTR1180 

TRTR1190 

TRTR1200 

TRTR1210 

TRTR122G 

TRTR123C 

TRTR124C 

TRTR1250 

TRTR1260 

TRTR1270 

TRTR128C 

TRTR129C 

TRTR1300 

TRTR1310 

TRTR1320 

TRTR1330 

TRTR1340 

TRTR1350 

TRTR136D 

TRTR1370 

TRTR138C 

TRTR1390 

TRTR14DC 

TRTR1410 

TRTR142C 

TRTR143C 

TRTR144C 

TRTR145C 

TRTR146D 

TRTR1470 

TRTR148D 

TRTR149G 
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15C 

C 

TRTR15CC 

151 

CALL OUPOUTC-iilCGf JSAVEf USD 

TRTR1510 

152 

RETURN 

TRTR152C 

153 

IICO PRINT llGlf NSS 

TRTR1530 

154 

1101 FORMAT! • SUBROUTINE CALLED FROM TRINTR. NS$-*I4) 

TRTR154C 

155 

RETURN 1 

TRTR1550 

156 

END 

TRTR156G 
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TREA»FTRE£.PATHIN 


1 

C 


SUBRCUTINO PATHIN ZEROS OUT ALL NUMBER BLOCKS IN 

ISAVE 

THAT ARE 

PAIN 

1C 

2 

C 


NOT EQUAL TO ANY OF THE N3IN ARRAY NUMBER BLOCKS. 

IF PATH HAS 50 

PAIN 

20 

3 

C 


ZEROES# STORE PATH NUMBER IN NPUP AND CALL PATHUP 

TO ELIMINATE. 

PAIN 

3C 

4 

C 


WRITTEN 3Y RON ECKSTEIN 1/1/70 



PAIN 

40 

5 

C 





PAIN 

50 

6 



SUBROUTINE PA THIN < $ # IS AV £ # JSl > 



PAIN 

60 

7 

C 





PAIN 

70 

3 



COMMON/ALLSUb/Mlf M2#H3fM 4# M 5# M6# M7 » LI #L 3#L 4 #L5» L6 



PAIN 

80 

9 



C0MM0N/SM5UB/NPUF(200 )#NPUPMX»NeCUT(5C) » NBOT MX # NBIN ( 50 1 

vNBINMXt 

PAIN 

90 

1C 



1 NSNUM(SG) fNBNMAX 



PAIN 

100 

11 



DIMENSION ISAVEIMl# M2) 



PAIN 

lie 

12 

C 





PAIN 

120 

13 

C 


INITIALIZE 



PAIN 

13C 

14 

C 





PAIN 

14G 

15 



DO 10 I-1*M2 



PAIN 

150 

16 



NPUP( I)=D 



PAIN 

160 

17 


10 

CONTINUE 



PAIN 

17C 

18 



L=C 



PAIN 

180 

19 

c 





PAIN 

190 

20 



00 IOC I'lfJSl 



PAIN 

200 

21 



ICOUNT=0 



PAIN 

210 

22 



DC 90 J^ltHl 



PAIN 

220 

23 



DO 80 K=1#NBINMX 



PAIN 

230 

24 



IF( ISAVEC J f I) .EQ.C) GO TO 85 



PAIN 

240 

25 



IF (NBIN(K).Eti. ABS( ISAVE (J #1) ))60 TO 90 



PAIN 

250 

26 


30 

CONTINUE 



PAIN 

260 

27 



ISAVE( JfDrO 



PAIN 

270 

28 


85 

IC0UNT=IC0UNT+1 



PAIN 

280 

29 


90 

CONTINUE 



PAIN 

290 

30 

c 





PAIN 

300 

31 

c 


IF 50 ZEROES STORE PATH NUMBER IN NPUP. 



PAIN 

31C 

32 

c 





PAIN 

320 

33 



IF (IC0UNT.lt. M l) 60 TO ICC 



PAIN 

330 

34 



L-L+1 



PAIN 

340 

35 



NPUP(L)-I 



PAIN 

350 

36 


100 

CONTINUE 



PAIN 

360 

37 



NPUPMX=L 



PAIN 

370 

38 



IF( NPUPMX.EQ.G) GO TO 1000 



PAIN 

3 80 

39 



CALL PATHUP($110C»ISAVE #JS1) 



PAIN 

39C 

40 

1000 

RETURN 



PAIN 

400 

41 

1100 

PRINT llCl 



PAIN 

41C 

42 

1101 

FORMATi* ERROR IN PATHUP CALLED FROM PATHIN.*) 



PAIN 

420 

43 



RETURN 1 



PAIN 

430 

44 



END 



PAIN 

440 
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trea»ftree;.trprnt 


1 

C 


subroutine; TRPRNT prints tree PATHSt ARRAY! ANO LOCATION OF 

CALL. TRPT 

10 

2 

C 


WRITTEN BY RON ECKSTEIN 1/1/70 

7RPT 

20 

3 

C 



TRPT 

3G 

4 



SUBROUTINE T RPRNT ( 1 1 ISA VE tJSltNSEQ ) 

TRPT 

40 

5 

c 



TRPT 

50 

& 



C0MM0N/ALLSUB/M1«M2»M5#M4 t M5 f M 6 f M 7 f L 1 1 L3 ? L4 * L5 t L6 

TRPT 

60 

7 



OIMENSiON ISAVECMl* M2) 

TRPT 

70 

6 

c 



TRPT 

80 

9 



GO TO (lf2#3f 4f2u) 

TRPT 

90 

10 


1 

PRINT 5» NSEQ 

TRPT 

IOC 

11 


5 

FORMAK •OISAVE AT ’IS) 

TRPT 

110 

12 



GO TO 2D 

TRPT 

120 

13 


2 

PRINT N5E3 

TRPT 

130 

14 


6 

FORMAK *CJSAVE AT ’IS) 

TRPT 

14C 

15 



GO TO 20 

TRPT 

150 

Id 


3 

PRINT 1 * NSEG 

TRPT 

160 

17 


7 

FCRMATJ *0K5AV£ AT *15) 

TRPT 

170 

18 



GO TO 20 

TRPT 

180 

19 


4 

PRINT 3» NSEO 

TRPT 

190 

20 


8 

FORMAT* ’CLSAVE AT ’IS) 

TRPT 

20 0 

21 


20 

DO 100 J-1»JS1 

TRPT 

210 

22 



DO 120 N=itMi 

TRPT 

220 

23 



IF* ISAVEt ^4# J) .EQ.O) GO TO 130 

TRPT 

221 

24 


120 

CONTINUE 

TRPT 

222 

25 



NzMl+1 

TRPT 

223 

26 


13C 

lEND^N'l 

TRPT 

224 

27 



PRINT 110# (ISAV£(N# J) #N:: ItlENO) 

TRPT 

23C 

28 


110 

FORMAT ( 1HD25 iI3 fix) *3(/lH 25 (I3#1XJ) ) 

TRPT 

24C 

29 


IGC 

CONTINUE 

TRPT 

2 50 

3C 



RETURN 

TRPT 

260 

31 



END 

TRPT 

2 70 


26 
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TREA»FTRE£.TRE£ 


1 

C 

SUBROUTINE TREE GENERATES THE PROBABILITY TREE SUCCESS PATHS. 

TREE 

1C 

2 

C 

WRITTEN BY PAUL CHELSON AND REVISED BY RON ECKSTEIN 1/1/70 

TREE 

20 

3 

C 


TREE 

3C 

4 

C 

SUBROUTINE TREE INPUTS IN COMMON. 

TREE 

40 

5 

C 


TREE 

5C 

6 

C 

IB(15»2»50) = IBdflfBLOCK NUMBERi^MAXlMUM NUMBER OF INPUTS OF NB 

.TREE 

60 

7 

C 

IB(lf2tBLGCK NUMBER >=MAXIMUM NUMBER OF OUTPUTS OF NBTREE 

70 

8 

C 

THE REST OF THE ARRAY = INPUTS AND OUTPUTS OF BLOCK NUMBERSCNB) 

TREE 

80 

9 

C 

ISUC = SUCCESS BLOCK NUMEER(BASt OF PROBABILITY TREE). 

TREE 

90 

10 

c 

JSUC = ISUC TO PRINT ORIGINAL PROS TREEt THEN ISUC VARIES. 

TREE 

ICO 

11 

c 

IPRINT = PRIM ALL TREES IF = It PRINT ORIGINAL TREE ONLY IF = 0. 

TREE 

lie 

12 

c 

NOPRNT z 1/C = OO/DO NOT PRINT TREES WHEN MODIFYING PROS DATA. 

TREE 

120 

13 

c 

IDR(5C) z HOLDS BLOCKS MADE INACTIVE BY A FAILED BLOCKIID). 

TREE 

13 C 

14 

c 

N1 z MAXIMUM NUMBER OF EQUIVALENT BLOCK SETS FROM MAIN PROGRAM. 

TREE 

140 

15 

c 

ITEMP z ARRAY OF EQUIVALENT BLOCKS FROM MAIN PROGRAM. 

TREE 

150 

1& 

c 


TREE 

160 

17 

c 

SUBROUTINE TREE OUTPUTS IN COMMON. 

TREE 

170 

18 

c 


TREE 

180 

19 

c 

ISAVE(5Ct2C0 ) - SUCCESS PATHS t + - BLOCK NUMBERt PATH NUMBER) 

TREE 

19 0 

20 

c 

J51 = MAXIMUM NUMBER OF SUCCESS PATHS. 

TREE 

200 

21 

c 


TREE 

210 

22 

c 

subroutine TREE VARIABLES. 

TREE 

220 

23 

c 


TREE 

23C 

24 

c 

IS = BLOCK NUMBER BEING WORKED ON. 

TREE 

240 

25 

c 

JS z path number being WORKED ON. 

TREE 

250 

26 

c 

ID - A FAILED BLOCK. 

TREE 

2 60 

27 

c 

lEl r NUMBER OF INPUTS FOR EACH OUTPUT BLOCK. 

TREE 

27 C 

28 

c 


TREE 

2 80 

29 


SUBROUTINE TREE(i) 

TREE 

29 0 

30 

c 


TREE 

300 

31 


CCMM0N/ALLSU&/KliM2fM3fM4iM5fM6tM7fLltL3fL4tL5tL6 

TREE 

310 

32 


C0MM0N/LGSUB/IB(15f 2f 50) t IIS ( 50) t IRB < 50 1 30 J • IT£MP< 4 1* 20 ) # N1 # 

TREE 

3 20 

33 


1 ISAVE(5Ct20Q) t ISUC t JSU C t JSl t JS AVE ( 50 t 20 C ) t I JSl t TOT t PROS « 65 ) • 

TREE 

330 

34 


2 PR0ai(65) tlPRINT.NSTOBY (15) t NSBYMXt NOPRNT 

TREE 

340 

35 


DIMENSION IDk(50) 

TREE 

35 C 

36 


DOUBLE PRECISION TO Tt PRO Bt PRO 81 

TREE 

360 

37 

c 


TREE 

37C 

38 

c 

INITIALIZE 

TREE 

380 

39 

c 


TREE 

39 0 

40 


JDZO 

TREE 

400 

41 


DO 2D IzitMl 

TREE 

41C 

42 


IDRtI)=0 

TREE 

420 

43 


DO 2C Kz1,M2 

TREE 

430 

44 

20 

ISAVE(IfK)=G 

TREE 

440 

45 


ITRYzlSUC 

TREE 

450 

46 


JSzl 

TREE 

460 

47 


IS = 1 

TREE 

470 

48 

c 


TREE 

4 60 

49 

.700 

ISAVEtlSf JS)zlTRY 

TREE 

490 

50 

701 

KzISAV£( IStJS ) 

TREE 

500 

51 


IF (K)7C3t7D2 i7D3 

TREE 

510 

52 

703 

IS=IS+1 

TREE 

520 

53 

C 


TREE 

530 

54 

C 

CHECK TO SEE IF IN DEAD ARRAY 

TREE 

54C 

55 

C 


TREE 

550 

56 


jzl 

TREE 

560 

57 

3 

JZj+1 

TREE 

570 

58 


DO 1 I=ltMl 

TREE 

5 00 

59 


IF(IE(JfltK)-'lDR(I))lt2tl 

TREE 

590 

60 

2 

IF(I3< Jfl»K))3t4f 3 

TREE 

600 

61 

1 

CONTINUE 

TREE 

510 

62 

4 

ISAVE(IStJS)zIB(Jtl#K) 

TREE 

620 

63 


GO TO 701 

TREE 

6 30 

64 

702 

ISl=IS-2 

TREE 

64 0 

65 

C 


TREE 

650 

66 

c 

CHECK TO SEE IF MAXIMUM NUMBER OF PATHS REACHED. 

TREE 

660 

67 

c 


TREE 

670 

68 


IF(JS“ M2iy02'»9Clt9D2 

TREE 

680 

69 

901 

IDR1=9C2 

TREE 

630 

70 


PRINT eODt IDRl 

TREE 

700 

71 

900 

FCRMAT< dERROR IN SUBROUTINE TREE AT*I4) 

TREE 

710 

72 


PRINT 910* JS 

TREE 

72 0 

73 

910 

FORMATdCJSt THE NUMBER OF SUCCESS PATHS GENERATED* EQUALS M2* THETREE 

730 

74 

1 MAXIMUM NUMClR OF SUCCLSS PATHS. •/• THE DIAGRAM SHOULD BE BROKEN 

TREE 

74 C 


75 20CUN INTO SEVERAL SUS-DIAGRAMS CONNECTED IN SERIES OR PARALLEL ANDTREE 750 
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76 


RERUN INDIVIDUALLY.* /• THEN ENTER 

EACH 

SUB-DIAGRAM AS ONE 

BLOCK 

TREE 760 

77 

40N A FINAL RUN.*/* 00 NOT BREAK 

UP 

EQUIVALENT 

BLOCK SETS OR 

STAN08TREE 770 

78 

5Y REDUNDANT SYSTEMS. */* IF Y CUR 

D1A6RAH CANNOT 

BE 

BROKEN DOWN SUFFTREE 7SC 

79 

SiCIENTLYf CONTACT PAUL CHELSON OR RON ECKSTEIN 

AT 

213 354-5691.*/ 

•TREE 790 

80 


7 JS=*I4> 







TREE 80C 

81 


RETURN 1 







TREE 810 

82 

C 








TREE 820 

83 

C 

CHECK PATH TO SEE IF ELIMINATED 

BY 

EQUIVALENT 

BLOCKS IN IT. 


TREE 830 

84 

C 

CIF THE SIGN OF ANY TWO EQUIVALENT 

BLOCKS 

IN A 

SET 

ALTERNATES! 

TREE 840 

85 

C 








TREE 850 

86 

90 2 

IF<N1.EQ.C» GO TO 1712 







TREE 860 

87 


00 1150 L-lfM7 







TREE 870 

88 


ITEMP(1»LU0 







TREE 68C 

39 

1130 

CONTINUE 







TREE 8 90 

90 


DO 1711 J::itMl 







TREE 9CC 

91 


K-ISAVE< J» JS) 







TREE 910 

92 


IFCKlll08tl7lifil08 







TREE 920 

93 

1108 

IF(Nl)l7llfl711»ll07 







TREE 930 

94 

1107 

DC HDD L-lfM 







TREE 940 

95 


DO 1100 M-2fL6 







TREE 950 

96 


IF CK-ITEMP(MtL> JllCCtllCl til CO 







TREE 960 

97 

HOC 

CONTINUE 







TREE 970 

98 


GC TO 1711 







TREE 980 

99 

1101 

IF( ITEMPdtLJ )lil2tll0 3t 1104 







TREE 990 

IDG 

1103 

ITEMPClfLUK 







TREEIOOO 

101 


GO TO 1711 







TREElOlO 

1C2 

1104 

IF (K>1102tl7iXtl711 







TREE1C2C 

103 

1112 

IF( K)l7lltllC2#1102 







TREE1030 

104 

1102 

ISAVEdS-lt JS>--ISAVECIS-'ltJSJ 







TREE1C40 

105 


JS=JS-1 







TREE1C50 

106 


GO TO 51G 







TREE1060 

107 

1711 

CONTINUE 







TREE1070 

108 

C 








TREE1C80 

109 

1712 

DO 500 N-lfISi 







TREE109C 

110 

500 

ISAVEiNt JS<d UISAVE(NfJS) 







TREEllOO 

111 


ISAV£( IS-it JS+1)--ISAV£( IS-lt JS) 







TREElllO 

112 

C 








TREE112C 

113 

c 

IO=A FAILED oLOCKtlOR HOLDS FAILED 

BLOCKS 

AND 

BLOCKS MADE 


TREE1130 

114 

c 

INACTIVE BY A FAILED SLOCK. AN INACTIVE BLOCK 

IS A 

BLOCK FOR 

UHICHTREE1131 

115 

c 

ITS INPUTS ARE ALL FAlLEDt THUS 

IT 

IS OF 

NO USE AS 

A POSSIBLE 

TREE1132 

116 

c 

BLOCK IN A SUCCESS PATH 







TREE1133 

117 

c 








TREE1140 

118 

510 

IO=-ISAVE( IS-lf JS^ll 







TREE1150 

119 


JS=JS-H 







TREE1160 

120 


JD3=0 







TREE1170 

121 

C 








TREE1180 

122 

C 

BEGIN DEAD ROUTINE WITH ID 







TREE1190 

123 

C 

CHECK ALL BLOCKS FOR WHICH ID IS 

AN 

INPUT 

BLOCK. WHENEVER ALL 

TREE1191 

124 

C 

INPUTS TO A BLOCK ARE IN lORCTHUS FAILED 

OR INACTIVE! t THEN 

THAT 

TREE1192 

125 

c 

BLOCK IS ALSO iNACTIVEtSO IS PUT 

IN 

lOR. 





TREE1I93 

126 

c 








TREE12CC 

127 

397 

IDR( JOtlHID 







TREE1210 

128 


J01=l+J03 







TREE1220 

129 

380 

JD = JD-tl 







TREE1230 

130 


J01=JD1-1 







TREEI240 

131 


IF (IDR( J0n379t378t379 







TREE125C 

132 

379 

I0R1=IDR( JOI 







TREE1260 

133 


K=:I6(It2tIDRll 







TREE1270 

134 

C 








TREE1280 

135 

C 

K IS NO. OF OUTPUTS OF IDRtJDI 







TREE129C 

136 

C 








TREE13Q0 

137 


DO 35C l=lfK 







TREE1310 

138 


I0R2-IBC I+lt2 tIORl) 







TREE1320 

139 


IB1=IB( lfltICR2 ) 







TREE1330 

140 

c 








TREE1340 

141 

c 

IBl IS NO. OF INPUTS FOR EACH OUTPUT BLOCK 




TREE1350 

14 2 

c 








TREE1360 

143 


IFfIBl-2l352t353*353 







TREE137C 

144 

353 

DO 354 N^ltlBl 







TREE1380 

145 


JN=1 







TREE139C 

146 


N2=N*1 







TREE14C0 

147 


DO 355 MritKi 







TREE141C 

148 


IF( IBCN2tlf I0R2)-I0R(MJ) 355t356t 

355 






TREE142C 

149 

356 

JN-2 







TREE143C 

150 

355 

CONTINUE 







TREE1440 

151 


GO TO C35Ct354ltJN 







TREE1450 
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152 


354 

CONTINUE 




TREE146Q 

153 

c 






TREE1461 

154 

C 


FINISHING THE 354 LOOP MEANS 

THAT ALL INPUTS TO 

IDR2 ARE IN 

IDR. TREE1462 

155 

c 


THUS IDR2 SHOULD ALSO BE INACTIVE AND IN IDR 



TREE1463 

156 

c 






TREE1464 

157 

c 


IS IT ALREADY IN IDR 




TREE1466 

153 

Q 






TREE1467 

159 


352 

DC 390 N=ifMl 




TREE147C 

160 



IF( I0R2-IDRIN ) )39Cf 35Cf330 




TREE148C 

161 


390 

CONTINUE 




TREE149C 

162 

c 






TREE1491 

163 

c 


IF NOT» PUT IDR2 IN IDR 




TREE1492 

164 

c 






TREE1493 

165 



JC1= JDl+1 




TREE150C 

166 



JD2-J0* JOl 




TREE1510 

167 



IDR( JD2i>IDR2 




TREE1520 

168 



IF« lORl JD2)-ISUC)350f 381t350 




TREE153Q 

169 


350 

CONTINUE 




TREE1540 

170 



GO TO 330 




TREE155C 

171 

c 






TREE156C 

172 

c 


381 ROUTE INDICATES TOTAL FAILURE 



TREE1570 

173 

c 






TREE158C 

174 


381 

K-0 




TREE1590 

175 


395 

K=K^1 




TREE1600 

176 



IKl=Li-K 




TREE1610 

177 



IF(ISAVEtlKl.JS) J395t395t396 




TREE162C 

178 


396 

ISAVEdKl# JS)--ISAVE(IKI.JS) 




TREE163C 

179 



ID = --ISAVEClKlt JS) 




TREE164C 

180 



IK2-IK1+1 




TREE1650 

181 



IS:IK2 




TREE1660 

182 



DO 377 NrlK2#Kl 




TREE1670 

183 


377 

ISAVE IN»JS>=L 




TREE168C 

184 



IF< IK 1-1)3 98# 300 1 398 




TREE1690 

185 


398 

DO 361 I=ltMl 




TREE17CC 

186 


361 

IOR{I)=0 




TREE1710 

187 



JO:C 




TREE172C 

188 



IK3=IK1-1 




TREE173C 

189 



JD3=1 




TREE1740 

19D 



DO 611 I::lflK3 




TREE1750 

191 



IF(ISAVECI*JS) )612f613t611 




TREE176C 

192 

c 






TREE1770 

193 

c 


ERROR 613t A SUCCESS PATH HAS 

A ZERO INSIDE IT. 



TREE178C 

194 

c 






TREE179D 

195 


613 

ISUC’613 




TREEieCC 

196 



PRINT 90Q» ISUC 




TREE1810 

197 



PRINT 950 




TREE182C 

198 


950 

FORMAT! ’OCHECK THE DIAGRAM ON 

THE FIRST PAGE OF 

THE PRINT OUT TO STREE1830 

199 


lEE IF ALL BLOCK NUMBERS ARE LISTED AND THAT THE 

INPUT/OUTPUT 

LISTSTREE184C 

2C0 



2*/» CORRECTLY DESCRIBE YOUR DIAGRAM. IF SO# MAKE 

SURE YOUR 

DIAGRATREE185Q 

201 


3M SATISFIES THE DEFINITION OF 

A RELIABILITY BLOCK 

DIAGRAM.* ) 

TREE1860 

202 



PRINT 32Gf JS 




TREE187G 

2C3 


920 

FORMAT! •GA SUCCESS PATH HAS A 

ZERO INSIDE XT. 

JS; 

= *I4) 

TREE188C 

204 



RETURN 1 




TREE1890 

205 

c 






TREE19CC 

206 


612 

JD3-JD3+1 




TREE191Q 

207 



IDR! JD3 I--IS AVE! 1 1 JS J 




TREE192C 

208 


611 

CONTINUE 




TREE1930 

209 



JD3- JD3-1 




TREE134G 

210 



60 TO 337 




TREE1950 

211 


378 

J0= JD-1 




TREE1960 

212 



DO 385 I=l#Mi 




TREE1370 

213 



IirLl-I 




TREE1980 

214 



ISA=ISAVElIlfJS) 




TREE1990 

215 



IF! ISA)3S5#385#61C 




TREE20Q0 

216 


610 

DO 386 J=1»M1 




TREE2010 

217 



IF! ISA-IDR! J) ) 386 #385 #38 6 




TREE2020 

218 


386 

CONTINUE 




TREE2C3C 

219 



GO TO lOCQ 




TREE2040 

220 


385 

CONTINUE 




TREE2G5C 

221 

c 






TREE2C60 

222 

c 


ERROR 385i THE 365 DO LOOP HAS BEEN COMPLETED. 



TREE207C 

223 

c 






TREE2080 

224 



ISUC::385 




TREE2D9C 

225 



PRINT 9C0r ISUC 




TREE2100 

226 



PRINT 95D 




TREE211C 

227 



PRINT 930# JS 




TREE2120 
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228 


S3C 

F0RMAT(*0TH£ 385 DO LOOP HAS BEEN COMPLETED. 

JS=d4 1 



TREE213C 

229 



RETURN 1 





TREE2140 

230 

C 







TREE215D 

231 


iCOC 

JN=IB<lfl*ISA> 





TREE2160 

232 



DO 387 J=itJN 





TREE2170 

233 



N2=IB< J+lf IfISA) 





TREE2180 

234 



DO 368 





TREE2190 

235 



IF( N2-*I0R( I) ) 388 1 387 f 3 83 





TREE22GG 

236 


388 

CONTINUE 





TREE2210 

237 



GO TO 63C 





TREE2220 

238 


387 

CONTINUE 





TREE223C 

239 

C 







TREE2240 

240 

c 


ERROR 387 t THE 387 DO LOOP HAS 

BEEN COMPLETED 

. 



TREE225C 

241 

c 







TREE226C 

242 



ISUC=387 





TREE227C 

243 



PRINT 9Q0# ISUC 





TREE2230 

244 



PRINT 950 





TREE229C 

245 



PRINT 940# ISA 





TREE230C 

246 


940 

FORMAT! *GISAf BLOCK NUMBER N CW 

BEING PUT INTO 

PATH# IS NOT 

AN 

ELEMTREE231C 

247 



LENT OF lOR ARRAY# YET ALL ITS INPUTS ARE ELEMENTS. ISA=*I4) 


TREE2320 

248 



RETURN 1 





TREE2330 

249 


630 

ITRY=N2 





TREE2340 

250 

c 







TREE235C 

251 

c 


END DEAD ROUTINE. RETURN ITRY 





TREE2360 

252 

c 







TREE237C 

253 



GO TO 700 





TR££23aO 

254 

c 







TREE2390 

255 

c 


SCO ROUTE MEANS END OF TREE 





TREE2400 

256 

c 


ZERO OUT FIRST BLOCK NUMBER IN 

FAILURE PATH. 




TREE241C 

257 

c 







TREE242G 

258 


SCO 

ISAVEd.JSJ-L 





TREE243C 

259 



JS1=JS-1 





TREE2440 

260 

c 







TREE245C 

261 

c 


IF IPRINT .NE. 4 - 00 NOT PRINT 

TREES AGAIN WHEN MOOYFYING 

PROB 

TREE246C 

262 

c 


IF IPRINT = 2 - PRINT ORIGINAL 

TREE ONLY. 




TREE248C 

263 

c 


IF IPRINT = 3 - PRINT ALL TREES 

• 




TREE2490 

264 

c 







TREE25GC 

265 



IF< IPRINT. lE.1) go to ICGOO 





TREE2510 

266 



IF (ISUC.NE.JSUC) GO TO 2000 





TREE252C 

267 



GO TO 2GG1 





TREE2530 

268 

2C0C 

IFdPRlNT.LT. 31 GO TO ICOCD 





TREE254C 

269 

2001 

PRINT 633# ISUC 





TREE255C 

270 


633 

FORMAT! dPROcABILITY TREEd4 1 





TREE256C 

271 



CALL TRPRNT!5»ISAVE#JS1#0) 





TREE257C 

272 

iCGDC 

RETURN 





TREE2580 

273 



END 





TREE259C 


30 
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IREA »FTREE.SYSP 


1 

C 


SUBROUTINC SYSP CALCULATES THE SYSTEM PROBABILITY - RETURNS TOT. 

SYSP 

10 

2 

C 


WRITTEN 3Y PAUL CHELSON AND REVISED BY RON ECKSTEIN 1/1/70 

SYSP 

20 

3 

C 



SYSP 

30 

4 

C 


SUBROUTINE SYSP INPUTS IN COMMON. 

SYSP 

40 

5 

c 



SYSP 

50 

6 

c 


JSl - MAXIMUM NUMBER OF SUCCESS PATHS FROM TREE SUBROUTINE. 

SYSP 

60 

7 

c 


ISAVE(5Cf200) = (BLOCK NUMBERSt SUCCESS PATH NUMBER! FROM TREE. 

SYSP 

7G 

8 

c 


N1 - MAXIMUM NUMBER OF EQUIVALENT BLOCK SETS FROM MAIN PROGRAM. 

SYSP 

80 

3 

c 


ITEMP X ARRAY OF EQUIVALENT BLOCKS FROM MAIN PROGRAM. 

SYSP 

90 

IQ 

c 


PROBl :: PROBABILITY OF THE BLOCKS FROM MAIN PROGRAM. 

SYSP 

ICO 

11 

c 



SYSP 

lie 

12 

c 


SUBROUTINE SYSP VARIABLES. 

SYSP 

120 

13 

c 



SYSP 

130 

14 

c 


PROS = PROBl 

SYSP 

140 

15 

c 


K = ISAVECJ»I) - BLOCK NUMBER BEING WORKED ON. 

SYSP 

15C 

16 

c 


PTOT r PROBABILITY OF ONE SUCCESS PATHIPRODUCT OF PROBCKJ IN PATH)SYSP 

160 

17 

c 


TOT = THE SUM OF ALL PTCT ( SYSTEM PROBABILITY ). 

SYSP 

17G 

18 

c 



SYSP 

180 

19 



SUBROUTINE SYSP 

SYSP 

190 

20 

c 



SYSP 

200 

21 



CCMMON/ALLSUb/Hl tM2fM3»M4 tMS tMo #M 7 tL 1 1 L3 t L4 » L5 »L6 

SYSP 

21C 

22 



C0MM0N/L6SU8/IB(15#2*5C) fIIS( SO) »IRB(50»30) • IT£MP( 4 1. 20 ) • N1 » 

SYSP 

220 

23 



1 ISAVEC5G»200) f ISUCt JSUC# JSl . J SAVE ( 5 C » 20C ) # I JSl tT OT * PROB ( 65 ) # 

SYSP 

23C 

24 



2 PROBl (65) tIPRINTfNSTDBY (15)f NSBYMX fNOPRNT 

SYSP 

2 40 

25 



DOUBLE PRECISION PTO T # T OT t PR CB t PRO Bl 

SYSP 

25C 

26 

c 



SYSP 

260 

27 

c 


INITIALIZE. 

SYSP 

27C 

23 

c 



SYSP 

280 

23 



DC 122 I=:liL4 

SYSP 

29 C 

30 



PROB(I)-PROBKI) 

SYSP 

300 

31 


122 

CONTINUE 

SYSP 

31C 

32 



T0T=O.OO+O 

SYSP 

320 

33 

c 



SYSP 

33 C 

34 

c 


ADO PROBABILITIES OF PATHS. 

SYSP 

340 

35 

c 



SYSP 

35 C 

36 



DO 710 I-ltJSl 

SYSP 

360 

37 



PTOTn. 

SYSP 

370 

38 



DO 130 L=1»M7 

SYSP 

380 

39 



ITEMP(ltL)=0 

SYSP 

39 0 

40 


130 

CONTINUE 

SYSP 

400 

41 

c 



SYSP 

410 

42 

c 


MULTIPLY PROBABILITIES OF BLOCK NUMBERS IN A PATH. 

SYSP 

420 

43 

c 



SYSP 

430 

44 



00 711 J^lfMl 

SYSP 

440 

45 



K=ISAVE( Jtl> 

SYSP 

450 

46 



IF( K)lQ8f7lOtiC8 

SYSP 

460 

47 

c 



SYSP 

470 

48 

c 


CHECK PATH TO SEE IF ELIMINATED BY EQUIVALENT BLOCKS IN IT. 

SYSP 

4 80 

49 

c 


(IF THE SIGN OF ANY TWO EQUIVALENT BLOCKS IN A SET ALTERNATES) 

SYSP 

49 0 

50 

c 



SYSP 

5CG 

51 


1C8 

IFINl)10ofl06fl07 

SYSP 

510 

52 


107 

DO 100 L=1»N1 

SYSP 

520 

53 



DO 100 M=:2»LC 

SYSP 

530 

54 



IF( K-ITEMP (MtL) ) 10C#1G1» IOC 

SYSP 

540 

55 


100 

CONTINUE 

SYSP 

55 0 

56 



GO TO 106 

SYSP 

560 

57 


101 

IF (ITEMP(lfL) J102tl03flC4 

SYSP 

570 

58 


103 

ITEMP(1»L)=K 

SYSP 

5 80 

59 



60 TO 1C6 

SYSP 

530 

60 


104 

IFI K)105#lQ9f 1C9 

SYSP 

600 

61 


109 

PR0B(KI=1.0 

SYSP 

610 

62 



GO TO 712 

SYSP 

620 

63 


105 

K=-K 

SYSP 

630 

64 



PROBl K)-1.0 

SYSP 

640 

65 



60 TO 110 

SYSP 

650 

66 


102 

IFI K)lll»112t 112 

SYSP 

660 

67 


112 

PROB(K)=C.C 

SYSP 

670 

68 



GO TO 712 

SYSP 

680 

69 


111 

K=-K 

SYSP 

69 0 

70 



PROBl K)=0.0 

SYSP 

700 

71 



GO TO 110 

SYSP 

710 

72 


106 

IFlK)7l4#711f 712 

SYSP 

720 

73 


712 

PTOT=PTCT*PROBIK) 

SYSP 

730 

74 



GO TO 730 

SYSP 

740 

75 


714 

K=“K 

SYSP 

750 
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76 

110 

PT0Tr?T0T*<l.C-PR08(KI J 

SYSP 

7 60 

77 

730 

PROS (KJ^PPOBKK J 

SYSP 

770 

78 

711 

CONTINUE 

SYSP 

780 

79 

710 

TOmOT^PTOT 

SYSP 

790 

8G 


RETURN 

SYSP 

800 

81 


END 

SYSP 

810 


32 
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TR£A*FTR££.STBYPR 


1 

C 


FUNCTION STBYPR CALCULATES THE STANDBY TREE PROBABILITY FOR TIME 

SBPR 

10 

2 

c 


T1 THRU TTOT. 

S6PR 

20 

3 

c 


WRITTEN BY RCN ECKSTEIN 1/1/7D 

SBPR 

30 

4 

c 



SBPR 

40 

5 



FUNCTION STBYPR(ISAVE*JS1 tRQ t TL t R T t K T f JT i TTO T t NST EP f TLD ) 

SBPR 

50 

6 

c 



SBPR 

60 

7 



C0MMCN/ALLSUB/MliM2fM3iM4*M5 #M 6 tM 7 tL 1 » L3 » L4 t L5 t LG 

SBPR 

70 

8 



DIMENSION ISAVE(Mlf M2 ) # RO ( Ml ) » TL < Ml) t RT (L4 *101) # TL D( Ml ) 

SBPR 

30 

9 



DOUBLE PRECISION TOT t FT OT *DP RE C i S TEP 1 * S TEP2 

SBPR 

9C 

10 

c 



SBPR 

100 

11 

c 


INITIALIZE 

SBPR 

lie 

12 

c 



SBPR 

120 

13 



STEP1=-TTCT* (KT-1 J/NSTEP 

SBPR 

130 

14 



ST£P2::-TT0T*JT/NSTEP 

SBPR 

140 

15 



TCTrC. 

SBPR 

15C 

16 

c 



SBPR 

160 

17 



DO 100 I’lfJSl 

SBPR 

170 

18 



PT0T=1.0 

SBPR 

180 

19 



DO 90 J^liMl 

SBPR 

19 C 

20 



K-ISAV£( J#I) 

SBPR 

200 

21 



IF (K)7Dt90*8C 

SBPR 

210 

22 


70 

K=-K 

SBPR 

220 

23 



IF (K.6T*M1) GO TO £6 

SBPR 

230 

24 



DPREC=(ST£P1*TL(K) )+(ST£P2*(TLD«K)-TL<K)n 

SBPR 

240 

25 



PTOT^PTOT* (l.C-RC(K) •DEXP (DPREC ) ) 

SBPR 

25C 

26 



GO TO 90 

SBPR 

260 

27 

c 



SBPR 

27G 

28 

c 


CHECK FOR STANDBY BLOC K( 51- 65 ) / ( L1-L4 ) . 

SBPR 

280 

29 

c 



SBPR 

29C 

30 


. 80 

IF< K.GT.Mi) GO TO 85 

SBPR 

300 

31 



DPREC=CSTEPl*TL(K))4(STEP2^rTLD(KJ“TL(K) M 

SBPR 

310 

32 



PTOT=PTOT+RQi K)*DEXP(OPR£C) 

SBPR 

320 

33 



GO TO 9C 

SBPR 

330 

34 


85 

PTOTrPTOT»RT<K*JT) 

SBPR 

340 

35 



60 TO 90 

SBPR 

350 

36 


86 

PT0T=PT0T*(1.C-RT(K*JT)) 

SBPR 

360 

37 


9C 

CONTINUE 

SBPR 

370 

38 


100 

T0T=T0T+PT0T 

SBPR 

3 80 

39 



STBYPRrTCT 

SBPR 

39 C 

40 



RETURN 

SBPR 

400 

41 



END 

SBPR 

410 
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TREA<«FTR£E.OUPTR£ 


1 

C 

SUBROUTINE DUPTRE DUPLICATES ISAVE INTO JSAVE AND SETS IJS1=JS1. 

DP TR 

1C 

2 

c 

WRITTEN BY RON ECKSTEIN 1/1/70 

OPTR 

20 

3 

c 


DPTR 

3C 

4 


SUBROUTINE OU PTRE ( JSA VEi I JS It IS A V£ # USD 

OPTR 

40 

5 

c 


DPTR 

50 

6 


C0MM0N/ALL5UB/MlfM2#M3#M4f M5* M6f M7 # LI # L 3# L 4 t L5 f L8 

DPTR 

60 

7 


DIMENSION ISAVE(MltM2)f JSAVE tMlfK2) 

DPTR 

7C 

3 

c 


OPTR 

80 

9 

c 

DUPLICATE ISAVE INTO JSAVE. SET IJSl r JSl. 

DPTR 

90 

10 

c 


OPTR 

ICO 

11 


DO 1C I-lrM2 

DPTR 

110 

12 


DO ID J-lfMl 

OPTR 

120 

13 


JSAVEC JtI)-ISAVE( Jill 

DPTR 

130 

14 


10 CONTINUE 

DPTR 

140 

15 


I JS1= JSl 

DPTR 

15C 

16 


RETURN 

OPTR 

160 

17 


END 

DPTR 

17C 


34 
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TR£A»FTREt.STDBY2 


1. 

C 

SUBROUTINE ST0BY2 IS FOR 2 OR MORE STANDBY BLOCKS CONTROLLED BY A 

ST8Y 

10 

2 

C 

SENSE BLOCK. IT GENERATES THE STANDBY TREEt CALLS SUBROUTINE 

ST BY 

2G 

3 

C 

TRINTR WHICH PUTS THE STANDBY SL0CK<51-65) IN PLACE OF THE FIRST 

STBY 

30 

4 

C 

PATH OP THE STANDBY TREE IMBEDDED IN THE PATH OF THE ORIGINAL TREESTEY 

40 

5 

C 

THEN ELIMINATES ORIGINAL TREE PATHS THAT CONTAIN 2ND-JS1 STANDBY 

ST3Y 

50 

& 

c 

PATHS. THEN PUTS THE NEGATIVE STANDBY BLOCK ( 51-65 ) /(L1-L4 J IN 

ST BY 

6C 

7 

c 

AND ELIMINATES CCRRESPON DING PATHS. WRITE STANDBY TREE ON TAPE. 

STBY 

70 

8 

c 

WRITTEN BY RGN ECKSTEIN 1/1/70 

STBY 

80 

9 

c 


STBY 

90 

10 


SUBROUTINE STDBY2{$tNS) 

STBY 

10 0 

11 

c 


STAY 

110 

12 


COMMON/ALLSUE/Kl tM2 »M3t M4 tM5 tM6 t K 7 t L 1 » L3 t L4 t L 5 tL6 

STBY 

12C 

13 


C0MM0N/LGSU8/IB(15f 2f 50) t IIS( 50 J t IRB < 50 » 30 ) tlTEMPi 41f 20) t Ni t 

STBY 

130 

14 


1 ISAVE(5D t2CCi) f ISUC » JSU C » USl tJSAVE ( 50 1200 ) »I JSl t TCT t FROB C65 ) t 

STBY 

14 0 

15 


2 PROBK 55) #IPRINT#NSTOSY (15)# NSBYMXfNOPRNT 

STBY 

150 

16 


CCMM0N/SMSUB/NFUF(2CD)tNFUPMXfNBGUT(ED> f NBOTHX.N6IN(50> »NBINKX# 

STBY 

16 C 

17 


1 NBNUM(50) iNBNMAX 

STBY 

170 

18 


DOUBLE PRECISION TOT # PROB #PR C31 

STBY 

19 0 

19 

c 


STBY 

200 

2D 

c 

INITIALIZE 

ST BY 

210 

21 

c 


STBY 

220 

22 


DO 1C I=1#M1 

STBY 

230 

23 


NBIN( I)=0 

STBY 

240 

24 


10 CONTINUE 

STBY 

25 0 

25 


NSS=0 

STBY 

260 

26 

c 


STBY 

270 

27 

c 

GENERATE SENSE TREE 

STBY 

280 

28 

c 


STBY 

290 

23 


ISUC=NS 

STEY 

30 0 

30 


CALL TREE( $1100) 

STBY 

310 

31 

c 


STEY 

320 

32 

c 


STBY 

330 

33 

c 

GENERATE STANDBY TREE BY ZERCIN6 OUT ALL BLOCK NUMBERS IN ISAVE 

STBY 

34 0 

34 

c 

EXCEPT THOSE STANDBY BLOCK NUMBERS CONTROLLED BY THE SENSE BLOCK. 

STBY 

350 

35 

c 


STBY 

360 

36 

c 


STBY 

370 

37 


IEND=:IRB(NS#i)4l 

STBY 

38 Q 

38 


L=a 

STBY 

3 90 

39 


DO 6C I=2#IEND 

STBY 

400 

40 


L=L+1 

STBY 

410 

41 


NBIN(L)=IRe(NS»I) 

STEY 

420 

42 


60 CONTINUE 

STBY 

430 

43 


NEINMX=L 

STBY 

44 0 

44 

c 


STBY 

4 50 

45 

c 

SUBROUTINE PATHIN ZEROS OUT ALL NUMBER BLOCKS IN ISAVE THAT ARE 

STEY 

460 

46 

c 

NOT EQUAL TO ANY OF THE NBIN ARRAY NUMBER BLOCKS. 

STBY 

470 

47 

c 


STBY 

480 

48 


CALL PATHIN($11G0#I5AVE# JSl) 

STBY 

4 90 

49 

c 


STEY 

50 0 

50 

c 

SUBROUTINE SFLEFT SHIFTS PATHS IN ISAVE TO THE LEFT UNTIL NONZERO 

STBY 

510 

51 

c 

NB AT LEFT-MOST. FILL IN ZEROS AFTER SHIFTING. IF 50 ZEROS PUTS 

STBY 

52 0 

52 

c 

PATH NUMBER IN NPUP ARRAY AND CALLS PATHUP TO ELIMINATE. 

STBY 

5 30 

53 

c 


STEY 

54 0 

54 


CALL SFLEFT($llQCf ISAVE# JSl) 

STBY 

550 

55 

c 


STEY 

560 

56 

c 

SUBROUTINE OUPOUT STORES PATH NUMBERS THAT ARE DUPLICATES IN NPUP 

STBY 

5 70 

57 

c 

AND CALLS PATHUP TO ELIMINATE DUPLICATE PATHS. 

STEY 

58 0 

58 

c 


STBY 

5 90 

53 


CALL DUP0UT($11CC»ISAVE »JSi) 

STBY 

60 0 

SC 

c 


STBY 

610 

61 

c 

SUBROUTINE NEGOUT STORES FAILURE PATH NUMBERS ENDING WITH A 

STBY 

62 0 

62 

c 

NEGATIVE SIGN IN NPUP ARRAY# THEN CALLS PATHUP TO ELIMINATE. 

STBY 

630 

63 

c 


STBY 

640 

84 


CALL NEG0UT($11CQ#ISAV£# JSl) 

STBY 

650 

65 

c 


STBY 

66 0 

66 

c 

NOW WE HAVE STANDBY TREE IN ISAVE. 

STBY 

670 

67 

c 


ST BY 

680 

68 

c 

SUBROUTINE TRINTR PUTS THE STANDBY BL0CK(51-65) IN PLACE OF THE 

STBY 

6 90 

69 

c 

FIRST PATH OF THE STANDBY TREE IMBEDDED IN THE PATH OF THE 

STBY 

700 

70 

c 

ORIGINAL TREE IN JSAVE# THEN ELIMINATES ORIGINAL TREE PATHS THAT 

STBY 

710 

71 

c 

CONTAIN 2ND- JSl STANDBY PATHS. THEN PUTS THE NEGATIVE STANDBY 

STBY 

720 

72 

c 

BLOCKi 51-85)/ (L1-L4) IN AND ELIMINATES CORRESPONDING PATHS. 

STBY 

730 

73 

c 


STBY 

74 0 

74 


CALL TRINTR ($1100# JSAVE# IJSl# IS A V£ . JS 1 # NST D8Y #:nS BY MX# NSS) 

STBY 

750 

75 

c 


STBY 

76 0 
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76 

C 

SAVE STANDBY TREE AND ITS STANDBY SLOCK C 51- 65 )/ C L1-L4) • 

STBY 

770 

77 

C 


ST BY 

78C 

78 


WRITEilO) NSTOBYtNSBYMX) f JSl 

STBY 

790 

79 


DO SCO I=ltM2 

STBY 

80C 

80 


WRITEdO ) (ISAVEUiI) » It Ml) 

STBY 

810 

81 

90C 

CONTINUE 

STBY 

820 

82 


RETURN 

STBY 

8 30 

83 

1100 

PRINT 1101 

STBY 

84C 

84 

1101 

FORMAT(» SUBROUTINE CALLED FROM ST0BY2.M 

STBY 

850 

85 


RETURN 1 

ST BY 

860 

86 


END 

STBY 

870 
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